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This  research  considers  the  design  of  a  representative  end-to-end  guidance  and 
navigation  system  for  the  terminal  phase  rendezvous  of  two  spacecraft  in  the 
circular  restricted  three-body  problem.  It  eissumes  the  target  vehicle  is  in  a 
small  radius  translunar  halo  orbit;  the  target  vehicle  remains  passive  with  re¬ 
spect  to  the  rendezvous  profile.  The  chaser  vehicle  is  in  an  independent  small 
radius  translunar  halo  orbit  and  executes  ail  rendezvous  maneuvers.  First,  the 
halo  orbit  guidance  problem  is  formulated  in  the  frequency- domain  from  which 
an  output  feedback  guidance  law  is  developed  using  7^2  control  theory.  Linear 
simulation  results  validate  the  guidance  law  and  provide  data  which  quantify 
the  effect  of  control  inputs,  noise  characteristics  and  halo  orbit  characteristics 


vi 


on  the  steady  state  halo  orbit  station-keeping  costs.  Secondly,  a  two  space¬ 
craft  terminal  phase  rendezvous  targeting  law  is  derived;  nonlinear  simulation 
results  validate  the  targeting  law.  Using  these  results,  the  three-body  ter¬ 
minal  phase  rendezvous  problem  is  further  characterized  and  contrasted  with 
the  two-body  problem.  Thirdly,  a  rendezvous  navigation  filter  is  constructed 
to  supply  the  targeting  law  with  chaser  vehicle  state  information.  Nonlinear 
simulation  results  validate  the  filter  design  and  provide  data  to  assess  the  filter 
performance.  Lastly,  two  rendezvous  scenarios  are  examined  to  demonstrate 
the  complete  system  design. 
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Chapter  1 


Introduction 


“  ...for  the  next  century,  back  to  the  Moon,  back  to  the  future, 
and  this  time,  back  to  stay.  And  then  a  journey  into  tomorrow,  a 
journey  to  another  planet,  a  manned  mission  to  Mars.  ” 

President  George  Bush 

In  his  speech  on  July  20,  1989,  President  Bush  launched  what  has 
become  known  as  the  Space  Exploration  Initiative.  The  goals  outlined  in  his 
speech  were  to  return  man  to  the  Moon,  establish  a  perniaiieiil  preseiuc  on 
the  Moon  and  then  proceed  on  to  Mars.  The  concept  of  human  exploration 
of  Mars  did  not  originate  with  the  Space  Exploration  Initiative;  initial  mis¬ 
sion  scenarios  were  proposed  by  von  Braun  and  Ryan  [1]  in  1954.  The  Space 
Exploration  Initiative  just  rejuvenated  national  enthusiasm  for  a  manned  mis¬ 
sion  to  Mars.  Its  significance  was  further  strengthen  by  the  “Report  of  the 
Advisory  Committee  On  the  Future  of  the  U.S.  Space  Program”  [2].  One  of 
the  15  specific  recommendations  of  the  Advisory  Committee  was  to  establish 
the  long-term  goal  of  the  U.  S.  Space  Program  to  be  the  human  exploration 
of  Mars. 

In  support  of  this  project,  NASA  conducted  a  90-day  study  of  the 
main  elements  of  a  Human  Exploration  Initiative  (3).  Iiirliidwl  in  N.AS.A's 
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report  were  four  approaches  to  a  manned  mission  to  Mars.  Common  to  all 
four  approaches  was  the  deployment  of  a  telecommunications  satellite  in  halo 
orbit  about  the  translunar  equilibrium  point  to  support  farside  lunar  commu¬ 
nications.  Figure  1.1  shows  the  geometry  of  a  telecommunication  satellite  in 
translunar  halo  orbit. 

In  parallel  with  the  Advisory  Committee  and  at  the  request  of  the 
National  Space  Council,  the  Synthais  Group  weis  formed  to  specifically  exam¬ 
ine  the  Space  Exploration  Initiative.  Their  report  [4]  was  similar  to  the  NASA 
90-day  study  in  that  they  also  proposed  four  architectures  to  a  manned  mission 
to  Mars  depending  on  the  chosen  scientific  objectives.  Although  a  telecommu¬ 
nications  satellite  in  translunar  halo  orbit  was  not  specifically  mentioned  in 
any  of  their  four  architectures,  the  Synth«is  Group  did  reference  traveling  to 
and  from  the  cislunar  equilibrium  point  and  even  possible  staging  operations 
from  this  point.  Further,  they  validated  the  need  for  new  or  additional  ren¬ 
dezvous  techniques  if  a  equilibrium  point  was  to  be  used  for  staging  operations 
to  Mars. 

Navigation  capabilities  were  also  addressed  by  the  Synthesis  Group 
as  a  supporting  teclinology.  They  concluded  that  navigation  from  onboard 
systems  would  satisfy  the  real  time  navigation  requirements  during  criticid 
operations  such  as  rendezvous.  In  addition.  NASA's  lead  center  for  the  Space 
Exploration  Initiative  has  identified  translunar  equilibrium  point  navigation 
as  a  potential  study  issue  [5]. 

In  this  chapter,  the  basic  premise  of  this  research  is  pul  forth  and 
specific  goals  are  detailed.  Previous  work  ndaliiig  lo  this  research  is  then 
summarized  and  followed  by  an  brief  outline  of  this  maiiiiscripl. 
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Figure  1 .2:  System  Design  Continuum 


1.1  Research  Premise 


The  primary  objective  of  this  research  is  to  design  a  representative 
end-to-end  guidance  and  navigation  system  (see  Fig.  1.2)  for  the  terminal 
phase  rendezvous  of  two  spacecraft  in  the  three-body  problem.  Further,  this 
system  is  demonstrated  using  a  realistic  and  practical  application;  the  major 
characteristics  and  trends  associated  with  this  application  arc  detailed.  The 
application  chosen  was  the  rendezvous  of  two  spacecraft  in  translunar  halo 
orbit;  the  translunar  halo  orbit  was  chosen  because  it  combines  the  mission 
concepts  introduced  in  the  NASA  90-day  study  [3]  and  the  proposed  architec¬ 
tures  of  the  Synthesis  Group’s  report  (4).  Figure  1.3  illustrates  the  rendezvous 
scenario.  The  secondary  objective  of  this  research  is  to  develop  an  interac¬ 
tive  nonlinear  simulation  package  for  terminal  phase  rendezvous  guidance  and 
navigation  system  evaluation  in  the  three-body  problem. 

Specific  goals  have  been  established  for  this  research: 


1.  Develop  a  frequency- domain  approach  to  the  translunar  halo  orbit  guid¬ 
ance  problem. 


2.  Formulate  the  system  model  for  the  translunar  halo  orbit  guidance  prob¬ 
lem  in  a  general  manner  applicable  to  both  'H2  and  Hoo  control  theory 
design  techniques. 

3.  Investigate  the  effect  of  system  parameters  on  the  resultant  translunar 
halo  orbit  geometry  and  station-keeping  cost. 

4.  Develop  a  two  spacecraft  terminal  phase  rendezvous  targeting  law  for 
the  three-body  problem. 

5.  Characterize  the  three-body  terminal  phase  rendezvous  problem  and 
contrast  it  with  the  well-known  two-body  problem. 

6.  Develop  a  rendezvous  navigation  filter  to  supply  chaser  state  information 
to  the  rendezvous  targeting  law. 

7.  Characterize  navigational  aspects  of  the  three-body  terminal  phase  ren¬ 
dezvous  problem. 

8.  Develop  a  nonlinear  simulation  package  to  be  used  interactively. 

9.  Construct  the  simulation  package  in  a  modular  manner. 

10.  Fully  document  all  simulation  source  code. 

11.  Provide  simulation  processing  status  messages  at  regular  intervals. 

1.2  Previous  Work 

1.2.1  Halo  Orbit  Guidance 

In  the  late  1960’s,  NASA  began  examining  follow-on  lunar  explo¬ 
ration  opportunities  to  the  the  Apollo  program.  In  an  effort  to  open  up  the 


farside  of  the  Moon  to  exploration,  studies  were  undertaken  to  determine  if 
continuous  communications  or  staging  operations  with  a  farside  lunar  base 
could  be  accomplished.  Three  main  concepts  were  examined: 

1.  A  translunar  halo  orbit. 

A  translunar  halo  orbit  is  an  pseudo-circular  orbit  about  the  translunar 
equilibrium  point  which  has  continuous  line-of-sight  with  the  Earth.  The 
term  “halo”  comes  from  the  trajectory  as  seen  from  the  Earth;  the  orbital 
path  looks  like  a  halo  around  the  Moon  (see  Fig.  1.1 ). 

2.  A  hummingbird  orbit. 

A  hummingbird  orbit  is  a  stationary  orbit  at  an  offset  point  from  the 
translunar  equilibrium  point.  This  orbit  is  similar  to  a  hummingbird  in 
stationary  flight. 

3.  A  polar  lunar  orbit. 

This  trajectory  is  a  classical  polar  orbit  around  the  Moon. 

Many  investigators  [6]  -  [26]  have  examined  the  general  halo  orbit 
guidance  problem;  this  section  will  summarize  only  those  sources  which  di¬ 
rectly  apply  to  the  translunar  halo  orbit  guidance  problem. 

Farquhar  [12]  [13]  and  Farqiihar  and  Kamel  [1.5]  provided  an  exten 
sive  review  of  previous  work  done  on  controlling  an  orbit  about  an  equilib¬ 
rium  point.  They  showed  that  using  a  simple  proportional  plus  derivative 
controller  provided  asymptotic  stability  while  minimizing  the  control  acceler¬ 
ation  required.  Deviations  from  the  desired  orbit  were  not  considered  in  the 
minimization.  Breakwell,  et  al  [7]  formulated  the  halo  orbit  guidance  problem 
as  a  periodic  system.  They  used  the  classical  optimal  control  approach  with 
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the  addition  of  an  observer  to  the  system  model.  Position  deviations  were 
considered  in  their  problem  formulation,  but  only  results  for  a  large  halo  orbit 
radius  were  given. 

Subsequently,  the  General  Electric  Company  completed  a  flight  dy¬ 
namics  study  [24]  of  both  the  halo  and  hummingbird  concepts.  They  found 
both  concepts  feasible,  but  preferred  the  halo  orbit  due  to  less  propulsion  re¬ 
quirements  for  station-keeping.  Their  formulation  used  a  frequency  matching 
guidance  law  with  discrete  impulses  applied  twice  an  orbit.  Heppenheimer  [18] 
used  phase-plane  methods  to  construct  a  family  of  locally  fuel-optimal  out- 
of-plane  period  controls.  More  recently,  Fraietta  and  Bond  [17]  computed 
station-keeping  costs  for  orbits  about  both  the  cislunar  and  translunar  equi¬ 
librium  points.  Vonbun  [26]  investigated  using  a  hummingbird  orbit  rather 
than  a  halo  orbit  and  found  in  general  it  required  ten  percent  more  accelera¬ 
tion  to  maintain  the  desired  position. 

Finally,  Farquhar  [14]  compared  the  use  of  a  polar  lunar  orbit  and 
halo  orbit  for  lunar  exploration  staging  operations.  He  concluded  that  a  halo 
orbit  space  station  could  offer  important  operational  and  performance  advan¬ 
tages  compared  to  a  polar  lunar  orbit  station.  Among  these  advantages  were 
increased  communication  opportunities  with  the  lunar  surface  and  inrreasv^d 
launch  windows  for  transfers  between  the  space  station  and  the  lunar  surface. 


► 
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1.2.2  Rendezvous  Targeting 

Rendezvous  with  a  spacecraft  in  translunar  halo  orbit  consists  of  a 
preliminary  and  a  terminal  phase.  The  preliminary  phase  guides  the  chaser 
vehicle  from  its  origin  to  the  halo  orbit  in  relatively  close  proximity  to  the 
target  vehicle.  A  typical  final  state  to  the  preliminary  phase  places  the  chaser 
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vehicle  15-20  km  behind  the  target  vehicle  with  its  velocity  vector  pointing  in 
the  direction  of  the  target  vehicle.  Subsequently,  the  terminal  phase  guides  the 
chaser  vehicle  to  a  small  specified  relative  position  and  velocity  with  respect 
to  the  target  vehicle.  The  Space  Shuttle  program  targets  terminal  phase 
rendezvous  to  place  the  chaser  vehicle  300  m  in  front  of  the  target  vehicle 
with  zero  relative  velocity. 

This  research  does  not  consider  the  preliminary  phase.  Many  inves¬ 
tigators  [27]  -  [34]  have  already  examined  how  to  optimally  guide  a  spacecraft 
from  low-Earth  or  low-lunar  orbit  to  the  translunar  equilibrium  point  or  to 
a  halo  orbit  about  the  translunar  equilibrium  point.  By  contrast,  published 
literature  on  terminal  phase  rendezvous  targeting  has  been  limited  solely  to 
the  two-body  problem. 

1.2.3  Rendezvous  Navigation 

The  rendezvous  targeting  law  is  dependent  on  accurate  state  infor¬ 
mation  of  the  chaser  and  target  vehicles.  Therefore,  a  rendezvous  naviga¬ 
tion  filter  is  required  to  provide  such  information.  .No  published  papers  w(‘re 
found  concerning  terminal  phase  rendezvous  navigation  in  the  three-body 
problem;  only  two  publications  were  found  on  the  orbit  determination  of  a 
single  spacecraft  about  an  equilibrium  point  in  the  Earth-Moon  three-body 
system  [24],  [35]. 

The  General  Electric  flight  dynamics  study  [24]  previously  cited  in 
the  halo  orbit  guidance  section  also  contained  an  initial  tracking  accuracy 
study  of  a  spacecraft  in  translunar  halo  orbit.  This  study  used  measurements 
from  Earth-based  sensors.  More  recently,  Gingiss  [35]  investigated  the  feasi¬ 
bility  of  using  current  radiometric  techniques  for  the  navigation  of  a  spacecraft 


to,  from,  or  at  the  Earth-Moon  equilibrium  points  using  various  navigation 
infrastructures,  but  not  including  onboard  navigation  systems.  These  two  pub¬ 
lications  showed  consistent  results  with  the  minor  differences  being  attributed 
to  the  advance  of  sensor  capabilities  in  the  last  20  years. 

1.3  Dissertation  Organization 

Cnapter  2  discusses  the  circular  restricted  problem  of  three- bodies. 
The  relevant  equations  of  motion  are  derived  and  the  location  and  stability 
of  the  equilibrium  points  determined.  Periodic  orbits  about  the  translunar 
equilibrium  point  are  discussed  and  the  need  for  active  halo  orbit  guidance 
established. 

Chapter  3  presents  the  halo  orbit  guidance  problem.  I'he  system 
models  are  developed  and  a  ^.2  optimal  guidance  law  computed.  A  linear 
simulation  is  formulated  to  validate  the  guidance  law.  Several  parametric 
studies  were  conducted  and  the  results  are  summarized. 

Chapter  4  develops  the  rendezvous  targeting  law.  The  target  and 
chaser  vehicle’s  state  equations  are  presented.  The  terminal  phase  initiation 
and  terminal  phase  final  targeting  laws  are  derived.  The  nonlinear  simulation 
described  in  Chapter  7  was  used  to  verify  this  law  and  the  results  from  several 
parametric  studies  are  given. 

Chapter  5  is  devoted  to  the  development  of  the  navigation  filter.  The 
dynamical  system  is  discussed  and  the  continuous-discrete  extended  Kalman 
filter  presented.  The  results  of  several  parametric  studies  are  given:  these  data 
were  generated  using  the  noiiliiiear  simulation  given  in  ('hapter  7. 

Chapter  6  contains  the  application  of  the  entire  guidance  and  nav¬ 
igation  system  to  the  two  spacecraft  rendezvous  problem  in  translunar  halo 
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orbit.  Two  specific  cases  are  examined;  the  first  case  is  an  entirely  “in-plane” 
relative  trajectory  for  the  rendezvous  profile  while  the  second  case  is  contains 
some  out-of-plane  relative  motion. 

Chapter  7  details  the  functional  description  of  the  nonlinear  sim¬ 
ulation  package.  This  simulation  package  was  developed  for  terminal  phase 
rendezvous  guidance  and  navigation  system  evaluation.  Each  module  associ¬ 
ated  with  a  mission  function  is  outlined. 

Finally,  Chapter  8  summaries  this  research  and  draws  general  conclu¬ 
sions  about  two  spacecraft  terminal  phase  rendezvous  in  the  circular  restricted 
three-body  problem. 


Chapter  2 


The  Circular  Restricted  Problem  of  Three-Bodies 


Investigation  into  the  problem  of  three-bodies  dates  back  to  the  eigh¬ 
teenth  century.  In  1772,  Lagrange  showed  the  existence  of  equilibrium  points 
in  the  restricted  three-body  problem.  The  three-body  problem  has  been  the 
focus  of  much  research  in  classical  celestial  mechanics;  .Szebehely  [.36]  devoted 
an  entire  volume  to  just  this  subject.  The  restricted  three-body  probieni  de¬ 
fines  a  system  where  two  primary  bodies  revolve  around  their  baryceiifer  in 
circular  orbits  under  the  influence  of  their  mutual  gravitational  attraction  and 
the  third  body,  which  has  significantly  less  mass,  is  attracted  by  the  two  pri¬ 
mary  bodies  but  not  does  not  influence  their  motion.  The  restricted  problem 
of  three-bodies  mathematically  describes  the  motion  of  the  third  body. 

Certain  approximations  are  made  when  using  the  restricted  problem. 
For  the  Earth-Moon  system,  Szebehely  (36)  makes  precise  statements  of  these 
approximations: 

1.  The  motion  of  the  primaries  (Earth  and  Moon)  takes  place  approxi¬ 
mately  in  a  plane; 

2.  The  motion  of  the  third  body  takes  place  i*sseiilially  in  the  same  plane; 

3.  The  effect  of  a  typical  probe  on  the  motion  of  the  primaries  is  of  tin' 
order  10"'®  &s  compared  to  the  forces  acting  between  the  primaries; 

4.  The  motion  of  the  primaries  is  essentially  circular. 
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The  first  approximation  is  satisfied  by  noting  that  the  oscillation  of 
the  Earth- Moon  orbital  plane  is  only  a  few  minutes  (1  min  =  0.0167  deg)  each 
year.  The  restricted  problem  used  in  this  research  will  be  modified  to  include 
three-dimensional  motion;  thus,  the  second  approximation  will  be  eliminated. 
The  third  approximation  restricts  the  third  body  to  less  than  13,000  lbs  which 
this  research  will  assume  to  be  satisfied.  Finally,  the  Moon’s  orbit  has  an 
eccentricity  of  0.055;  hence,  the  last  approximation  is  met. 

In  this  chapter,  the  nonlinear  equations  of  motion  are  developed 
in  both  dimensional  and  non-dimensional  units.  The  equilibrium  points  for 
the  Earth-Moon  system  are  identified  and  the  non-dimensional  equations  of 
motion  are  linearized  about  the  translunar  equilibrium  point.  Lastly,  the 
linearized  equations  are  solved  analytically  and  it  is  shown  that  a  small  radius 
halo  orbit  can  not  exist  without  active  guidance. 

2.1  Nonlinear  Equations  of  Motion 

Consider  the  motion  of  a  third  body  in  the  gravitational  field  of  two 
primary  bodies  of  spherically  symmetric  mass  distribution.  These  two  primary 
bodies  revolve  about  their  barycenter  in  circular  orbits.  Figure  2.1  shows  the 
force  diagram  for  this  problem.  .An  ineitial  coordinate  system  (i/.  }y.  kj)  and 
a  rotating  coordinate  system  (i/j,  j/j,  fc/j)  with  their  origins  at  the  oarycenter 
are  shown  in  Fig.  2.1.  The  relative  positions  of  the  third  body  with  respect 
to  the  Earth  and  Moon  are  given  by  Tgs  and  respectively.  Three  forces 
act  on  the  third  body:  the  gravitational  attraction  of  the  Moon  (Fa/),  the 
gravitational  attraction  of  the  Earth  (F/e),  and  the  control  force  generated  by 
internal  propulsion  (Fp).  The  positions  of  the  Moon  and  Earth  in  the  rotating 
coordinate  system  are  given  by  a  and  b.  The  orbital  motion  of  the  .Moon  is 
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Figure  2.1:  Force  Diagram 

described  by  the  mean  motion  (n)  and  dimensional  time  (<*). 

For  the  rotating  coordinate  system,  summing  the  external  forces 

yields 


=  Fm  +  Fe  +  Fp 


(GrmMmsTMs)  _  f CrmEmstES  1 

kMsI^  J  I  kESp  J 


kMsI^  j  I  kESp  J 

f  GrTTiMmslix'’  -  a)i/t  +  +  z’kfj]  1 

1  iTMsf  I 

f  GrmEmsKx'  +  b)»fl  -j-  ) 

1  J 

+  {mskxifi  +  Uyjit  + 


+  mscis 


where  (r*,y*,c’)  are  the  dimensional  po.sition  components  of  the  third  body 
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r 


A 

in  eaich  axis  of  the  rotating  coordinate  system  (*k,  Jr,  ib/i),  respectively.  Dif¬ 
ferentiating  rs  twice  in  the  rotating  coordinate  system  yields 

d^rs 


(2.4) 


Applying  Newton’s  second  law  of  motion  for  a  constant  mass  system 

(Prs 


=  ms 


dr^ 


(2.5) 


yields  the  three  second-order  scalar  differential  equations  of  motion  in  dimen¬ 
sional  coordinates 


o_“v  _2_.  ^  fmM(x’-a)  ,  mR(x*+b) 

~2n—-nx  =  p — ^  + 


<Px^  dy^ 
dr^  ^"^dt 


+  2n^-nV  =  -GrV 


df^ 


r  mM  .  rriE 


dt 


^  wm  ,  ms  1 


kfisl 

i|  +  ®y 


I  +  Or  (2. 


dt 


(2.6) 

(2.7) 

(2.8) 


where 


kiVfsl  =  \/(r'  -  a)*  +  y*^  + 
kesl  =  y/(r'  +  b)2  -|-  +  z‘^. 


(2.9) 

(2.10) 


The  non-dimensional  equations  of  motion  result  after  substituting 
the  non-dimensional  variables 


I  =  nt 


A  I 
X  = 


(a+b) 


(2.11) 

(2.12) 
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A 

y  = 


t/x 

u, 

t/. 


A 

A 

A 

A 

A 


y* 

(a+b) 


(a+b) 

rriM 

(m£  +  mAf) 
flr 

(a+b)n^ 


(a+b)n^ 

az 

(a+b)n^ 


(2.13) 

(2.14) 

(2.15) 

(2.16) 

(2.17) 

(2.18) 


into  Eqs.  (2.6  -  2.8) 


<Px 

w 

^dy  f/i[x-(l -/i)]  ,  (l-p)[x-|-p]'| 

dt  \  PMS^  pES^  j 

^  -h  t/x  (2.19) 

d^y 

de 

+  2^-,  =  -J  +  + 

dt  ( PMS  PES  ) 

(2.20) 

t  =  -4  "3+“‘S‘4+''. 

dt^  [PMS^  PES^  j 

(2.21) 

where 

PMS  =  -  (1 -/i)P +  y*  + 

(2.22) 

PES  =  \J[^  +  pY  +  y*  + 

(2.23) 

Equations  (2.19  -  2.21)  are  the  preferred  equations  of  motion  since  specific 
information  about  the  two  primary  bodies  (other  than  the  mass  ratio)  is  not 
required.  Also,  the  non-dimensional  equations  have  better  numerical  condi¬ 
tioning  since  all  terms  are  of  the  order  one. 


Define  the  nonlinear  state  and  control  as 

■  x(0  ■ 

z(0 

x(0 

y{0 
i(0 


^(0  = 


(2.24) 
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and 


Uii)  = 


LW 
Uyii) 

[  UAt)  J 


(2.25) 


Writing  Eqs.  (2.19  -  2.21)  in  first  order  functional  matrix  form  gives 

X  =  f,{X{t))+g,{U{t)) 


(2.26) 


where 


fiiXit))  = 


gduit))  = 


(•)  ^ 


y 

2 


X  + 


I  PMS  PES 

-  l] 

‘^XpMS  PES  J 


Pms 


3  + 


0 

0 

0 

Vr 

Uy 


di 


(2.27) 


(2.28) 


(2.29) 


2.2  Equilibrium  Point  Location 

Equilibrium  points  occur  when  all  external  forces  are  balanced.  Math¬ 
ematically,  they  occur  when  all  derivatives  in  the  differential  equations  are  zero 
or  undefined.  In  Eq.  (2.27),  setting  all  derivatives  to  zero  yields 

(l-p)(x-|-p) 


I  - 


Pms 


Pes 


=  0 


I  Pms  pes  J 


:  =  0. 


(2.30) 

(2.31) 

(2.32) 
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Equation  (2.31)  yields  two  cases  for  equilibrium  points; 


1 


1  - 


PMS^ 


y 

PES^  ) 


=  0 

=  0. 


(2.33) 

(2.34) 


The  first  case,  Eq.  (2.33),  locates  the  collinear  equilibrium  points;  the  second 
case,  Eq.  (2.34),  locates  the  triangular  equilibrium  points. 


Equation  (2.30)  was  derived  for  the  third  body  in  the  region  con¬ 
taining  the  Li  equilibrium  point.  For  the  third  body  in  the  regions  containing 
the  other  two  collinear  equilibrium  points,  the  signs  on  the  second  and  third 
terms  in  Eq.  (2.30)  change.  Thus,  Eq.  (2.30)  becomes 

/i(i-(l  -/i)|  (1  -fi)[x  +  n] 


XT 


=  0 


PMS 


PES 


(2.35) 


where  the  appropriate  signs  are  given  in  Table  2.1.  Equation  (2.35)  results  in 
a  fifth-order  equation  in  x  where  the  real  root  locates  the  equilibrium  point. 


Table  2.1:  Appropriate  Signs  for  Collinear  Equilibrium  Points 


Equilibrium  Point 

Second  Term 

Third  Term 

Lx 

- 

- 

L2 

- 

Lz 

+ 

+ 

For  the  triangular  equilibrium  points,  Eqs.  (2.30)  and  (2.34)  are 
satisfied  when 


PES  =  i  (•2.30) 

Pms  =  I-  (2.37) 

Hence,  Eqs.  (2.36)  and  (2.37)  are  used  instead  of  Eqs.  (2.30)  and  (2.34)  to 
locate  the  these  equilibrium  points. 
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For  the  Earth-Moon  system,  fi  =  0.0121506  [37],  which  results  in  the 
equilibrium  points  given  in  Table  2.2.  Figure  2.2  shows  the  geometry  of  the 
restricted  three-body  problem  equilibrium  points. 


Table  2.2:  Equilibrium  Point  Locations 


Equilibrium  Point 

X 

y 

z 

Translunar 

(^.) 

1.155682 

0 

0 

Cislunar 

(L2) 

0.836915 

0 

0 

Transearth 

(^3) 

-1.005063 

0 

0 

Triangular 

(f'4,5) 

0.487849 

±  0.866025  0 

2.3  Linearized  Equations  of  Motion 


and 


Define  the  linearized  state  and  control  as 


^  ^nom 

V 

y  ~~  J/nom 

c 

A 

^  ^nom 

i  ""  ^nom 

y  ~~  J/nom 

uJ 

.  ^  ^nom  , 

=  x-x 


nom 


A 

'  ' 

«  = 

— 

.  . 

u  -  n 

I  r  '  rnom 

I’  -  I' 

'■  V  Vnom 

V.  -  . 


U  -U 


tu>H4 


(2.38) 


(2.;19) 


where  (  denotes  a  nominal  value.  Neglecting  the  higher  order  terms  in 


a  Taylor  series  expansion  yields 


17  =  AqO  -|-  Bcu 


(2.40) 


where 


Ac 

Bg 


'?ll 

dx 


J  nom 


du 


(2.41) 


► 


t 


nom 


(2.42) 


Dj 


/  \ 

/  \ 

/  \ 


L2\  L, 


Earth 


Unstable 
L^.Lj  Stable 


Barycenter 


\  / 
\  / 

L, 


Moon 


Figure  2.2:  Fquilibriiim  Point  Locations 

Substituting  the  nominal  value  for  the  translunar  equilibrium  point  from 
Table  2.2  yields 


Ag  = 


0 

0 

0 

7.380861 

0 

0 


-2.190431 

0 


0  10  0 

0  0  10 

0  0  0  1 

0  0  2  0 

0  -2  0  0 
-3.190431  0  0  0 


(2.43) 


0  0  0 
0  0  0 
0  0  0 
I  0  0 
0  I  0 
0  0  I 


(•2.11) 
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Note  the  ^-axis  can  be  uncoupled  from  the  ^  and  7/-axes  after  the  linearization. 

2.4  Periodic  Orbits 

The  general  solution  to  Eq.  (2.40)  is 

=  4>g(<,0)  12(0)  +  /  -  7’,0)  Bg  u(r)  dr  (2.45) 

Jo 

where 

<Pg(<,0)  =  exp^°^  (2.46) 

and  satisfies 

<Pg(<,0)  =  >Ig<Pg(<,C).  (2.47) 

The  first  term  in  Eq.  (2.45)  represents  the  unforced  solution  to  Eq.  (2.40)  and 
the  second  term  gives  the  forced  solution.  Appendix  A.l  details  the  solution 
to  the  state  transition  m.atrix,  ^g(^O).  Hence,  the  unforced  solution  to  Eq. 
(2.40)  becomes 

^(<)  =  1.586694(:;,exp-''  -1.586694C'2exp^ 

-0.1 84326(73 cos(l. 8626470  -f  sin(  1.862647<)  (2.48) 

V{t)  =  C,exp-'>®«®^^‘+C2exp'‘®«"^'‘ 

+C4Cos(  1.8626470  +  7-r^7::sin(  1.8626470  (2.49) 

1 .862647 

C(0  =  C(0)cos(1.7861780  +  T-^^sin(1.7861780  (2.50) 

1.  tool  lo 

where 

Cl  =  0.420595^0) -0.12301^(0) 

+0.0786673r/(0)  +  0.077.5267//(0) 

Ci  =  -0.420595^(0) -0.12.301  U(0) 


(2.51) 
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-f0.07866737/(0)  -  0.077  52677/(0) 

(2.52) 

=  1.815860^(0)-!-  1.3347107/(0) 

(2..5.3) 

=  0.246022^(0) +  0.84  26657/(0). 

(2.54) 

In  general,  Eqs.  (2.48  -  2.50)  are  unbounded  because  of  the  exponentially 
increasing  terms;  thus,  the  linear  system  is  unstable.  However,  if  the  ini¬ 
tial  conditions  are  chosen  properly,  the  exponential  terms  can  be  eliminated. 
Selecting 

^(0)  =  0.6395147/(0)  (2.55) 

^(0)  =  -5.425170^(0)  (2.56) 

reduces  Eqs.  (2.48  -  2.50)  to  purely  periodic  terms 

^(0  =  ^(0)  cos(  1.862647/)  + 0.3433367/(0)  sin(  1.862647/)  (2.57) 

7/(/)  =  7/(0)cos(1.862647/)-  2.928719^(0)siii(1.862617/)  (2.58) 

C(()  =  C(0)cos(1.786178()  +  Y^^sin(l-7861780.  (2.59) 

In  addition,  if  the  periodic  orbit  is  chosen  to  be  centered  about  Z/j,  that  is 

^(0)  =  0,  and  the  initial  point  of  the  orbit  is  selected  so  that  ((0)  =  0, 
Eqs.  (2.57  -  2.59)  further  reduce  to 


i{t)  =  0.3433367/(0)  sin(1.862647/) 

(2.60) 

7/(/)  =  7/(0)  cos(  1.862647/) 

(2.61) 

C(/)  =  C(0)cos(  1.786178/). 

(2.62) 

It  is  readily  apparent  that  a  small  radius  halo  orbit  can  not  exist  without 
active  guidance  because  of  the  {liffereiit  natural  freqneiuies  betwi'en  the  // 
axis  and  the  (.'-axis.  Figure  2.3  graphically  illustrates  this  I’esnlt  by  clu)osing 


2;} 


7;(0)  and  ^(0)  to  be  0.0091 1  (3500  km)  and  90  deg  oul-of-phase.  I'liis  Lissajous 
trajectory  clearly  shows  periods  of  time  when  the  Moon  will  block  the  line-of- 
sight  between  the  Earth  and  a  third  body;  this  period  of  blockage  equates  to 
a  disk  with  a  radius  of  0.00806  (3100  km)  centered  at  the  origin  on  Fig.  2.3. 


^  (Non-dimensional) 
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Figure  2.3:  Uncontrolled  Periodic  Orbit 
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Chapter  3 


Halo  Orbit  Guidance  Law 


Several  classical  time-domain  techniques  have  been  investigated  for 
controlling  a  halo  orbit.  In  this  chapter,  a  frequency- domain  technique  will  be 
introduced.  "Wj  control  theory  applies  a  frequency-domain  formulation  to  the 
classical  time-domain  Linear  Quadratic  Gaussian  problem.  The  advantage  to 
the  frequency-domain  formulation  is  that  it  can  be  extended  to  modern  control 
theories,  such  as  'H2  and  Hoc  theory,  so  that  the  class  of  plant  disturbances 
and  measurement  noises  considered  can  be  expanded. 

This  chapter  starts  with  the  formulation  of  the  H2  problem  state¬ 
ment  and  system  models.  The  state  space  H2  solution  is  then  discussed,  A 
simulation  is  developed  to  validate  the  guidance  law  and  results  from  two 
typical  guidance  laws  are  presented.  Lastly,  results  from  several  parametric 
studies  are  given. 


3.1  Problem  Statement 


H2  control  theory  minimizes  the  2-norm  of  the  closed  loop  transfer 
function,  Tj^,,  between  the  regulated  variables.  Zre^,  and  the  exogenous  inputs, 
tr2!  the  H2  system  model  is  shown  in  Fig.  :LI.  I'liis  Iransfer  luncliun  i.s  the 
same  as  the  lower  linear  fractional  transformation  of  tin*  syst(*m  |)lant.  Pf,-. 
and  the  guidance  law,  Khaio-  Given 


Pg  = 


Pg\\  Pgu 
i’cji  P  an 


(3.1) 
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then 


Ttw  =  P GW  P G\2l^ kalo(I  —  P 022^^ halo)  ^  P G2\ 

=  F  t{P  G 1  K  halo)  ■ 

A  precise  problem  statement  is 
Minimize 

where  the  regulated  variables  are  chosen  to  be 


(3.2) 

(3.3) 

(3.4) 


^reg  — 


PnV 

PdC 

u 


(3.5) 


and  the  exogenous  inputs  are  chosen  to  be 


Wo  = 


Wd 


[3.6) 


The  plant  disturbances,  Wd,  and  the  measurement  noise.  tu,„.  have  a  fixed 
power  spectrum.  In  addition,  and  p(^  are  .selectable  constants  so  that  the 
relative  weighting  of  position  deviation  and  control  acceleration  can  be  varied. 
Since  the  deviation  in  the  ^-axis  does  not  affect  the  line-of-sight  between  the 
halo  orbit  and  Earth,  ^  was  not  included  as  a  regulated  variable. 

'H2  minimization  is  the  frequency-domain  equivalent  to  Linear  Quad¬ 
ratic  Gaussian  minimization  in  the  time  domain.  The  equivalent  problem 
statement  in  the  time-domain  is 


Minimize  ~  j  Q lqg^ ^ 


where 


Qiqg  -- 

9l.QGl  1  = 


0  0 

0  0  0 
0  0 
0  0  pcM 


(3.7) 


(3.N) 


(3.9) 


28 


► 


Plant  Disturbances  .  J 

(eg.  Nonlinear  Effects, 

Third  Primary  Body)  1 


Figure  3.2:  Desired  System  Model 


3.2  System  Models 


Figure  3.2  shows  the  desired  system  model  where  G  represents  the 
plant  and  Khaio  represents  the  guidance  law.  The  plant  and  guidance  law 
have  the  general  realizations 


and  halo  — 


Bk 

Cs 

Dk 

(3.10) 


1 


Figure  3.3  shows  an  expanded  system  model  where  G  has  l)eeii  l)roken  down 
into  its  realization  parts;  Gj-  has  the  realization 


G,  = 


-^G 


(3.11) 


For  this  research,  only  position  measurements  are  considered.  Hence, 

Cg  =  [  /  0  ]  and  Dg  =  0.  (3.12) 


I 


I 


The  vironment  consists  of  plant  disturbances  and  measurement  noise,  both 
mof  '  as  white  gaussian  noise  with  zero  mean  and  identity  fixed  power 


I 
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Figure  3.3:  Expanded  System  Model 
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spectrum.  The  plant  disturbance,  or  process  noise,  is  then  scaled  by  a  non- 
dimensional  matrix,  ilj. 


n, 


7r„/  0 

0  TT*/ 


(3.13) 


where  and  tt,  are  the  non-dimensional  square  root  of  the  fixed  power  spec¬ 
trum  in  velocity  and  acceleration,  respectively.  The  measurement  iioi.se  lias 
also  been  scaled  by  a  non-dimensional  matrix,  Um, 


TI  m  — 


(3.14) 


where  Xp  is  the  non-dimensional  square  root  of  the  fixed  power  spectrum  in 
position.  Figure  3.1  shows  the  system  model  required  by  Doyle,  et  al  [38]  for 
the  two  Ricatti  equation  "Hj  control  system  design  method.  The  plant  (Pc) 
in  this  system  model  has  the  general  realization 


A 

B, 

B2 

Pg  = 

c, 

Du 

D\2 

c. 

D21 

D22 

(3.15) 


and  requires  a  particular  structure  for  the  D  matrices.  Dailey  [39]  summa 
rized  a  method  for  loop  shifting  and  scaling  these  matrices:  this  method  uses 
the  singular  value  decomposition  of  the  Dij  and  D-n  malrict*s  to  transform 
a  general  system  model  into  the  desired  form.  Transforming  the  expanded 
system  model  (Fig.  3.3)  into  the  system  model  and  applying  scaling  yields 


(3.16) 


► 


I 


» 


» 


» 


;{i 

where  (  )  indicates  scaled  parameters  and  q,,  and  are  matrices  defined 
such  that 

t]  =  oirfO  and  ^  (3-17) 

3.3  State  Space  H2  Solution 

The  'H2  optimal  solution  involves  solving  two  algebraic  Ricatti  equa¬ 
tions  (see  Doyle,  et  al  [38]).  The  first  Ricatti  equation  (X2)  represents  state 
feedback;  the  second  Ricatti  equation  (Y 2)  represents  observer  feedback.  For 
this  system,  the  two  Ricatti  equations  are 

X2AG  —  X2BgBg^X2  +  Ac^ X2 Pc^oic^otc  —  ®  (3.18) 

Y2Ag'^ -Y2CG'^n„r'^n„r'CGY2  +  Ac.Y2-\-n,n/  =  o.  (3.19) 

Given  the  solutions  to  Eqs.  (3.18)  and  (3.19),  the  feedback  gain  matrix  be¬ 
comes 

Fni  =  -Bg^X2  (3.20) 

and  the  observer  gain  matrix  is 

(3.21) 

Finally,  the  optimum  guidance  law  is  given  by 


where 

Aii  =  A-  BgBg^X2  -  Y2CG‘'n,.-'^n^-'CG.  (3.24) 
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3.4  Simulation  Results 

A  simulation  was  developed  to  validate  the  guidance  law  given  in 
Eq.  (3.23).  For  the  initial  cases  examined,  linear  simulation  results  were 
found  to  be  in  agreement  with  nonlinear  simulation  results.  Therefore,  the 
linear  simulation  was  used  to  generate  the  results  presented  in  this  section 
and  Appendix  B.  The  state  equation  of  the  simulation  is  first  formulated; 
then,  data  from  two  typical  H-i  optimal  guidance  laws  are  given  followed  by 
data  from  three  parametric  studies. 

з. 4.1  State  Equation  Formulation 

The  mathematical  model  of  the  guidance  law,  Eq.  (3.23),  is 

Xk  =  +  (3.25) 

«  =  Cf^XK  (3.26) 

where  Xk  are  the  guidance  law  states.  The  inputs  to  the  guidance  law  are 
deviations  from  the  halo  orbit,  e,  and  the  outputs  are  the  control  accelerations, 

и.  From  the  expanded  system  model.  Fig.  3.3,  note 

e  =  Ccfi  +  +  r.  (3.27) 

Hence, 

Xk  =  +  A^Xk  +  +  B,-n,„  'r.  (3.2S) 

The  mathematical  model  associated  with  P  is 

h  = 


e 


AgB  +  n^Wd  +  Bc.il 

n„r^CaO  + 


(3.29) 

(3.30) 


Hence,  the  Kalman  filter  mathematical  model  yields 

h:=AGn  +  BGU-Ln2nm~'CGin-n)  (3.31) 

where  the  (  )  denotes  the  estimated  system  states. 

Assume  the  guidance  law  states  are  a  linear  combination  of  the  sys¬ 
tem  state  estimates  plus  additional  undetermined  states,  F.  Hence, 

x^  =  n  +  r.  (3.32) 

Differentiating  Eq.  (3.32);  substituting  Eqs.  (3.25),  (3.26),  (3.31)  and 
(3.32);  and  noting  u  =  u  yields 

r  =  {Ag  —  Bf^nm~^CG)r  —  +  Bf^nm~^T.  (3.33) 

The  final  state  equation  for  the  linear  simulation  is  formed  by  aug¬ 
menting  Eqs.  (3.28),  (3.29)  and  (3.33)  into  a  matrix  equation 


a 

Ag  BGCf; 

0 

n  ■ 

Xk 

r 

Bf^rijn  Co  Af^- 

0  0 

0 

{Ag-  Bf;n„r'CG) 

Xh 

r 

/7d  0  0 

+ 

0 

0  —L-H2  Bf^n„,~^ 

T 

and  noting 

a  =  XK-r  (3.35) 


3.4.2  Typical  'H2  Guidance  Laws 


The  simulation  was  run  for  300  days  (20  revolutions  of  the  halo  orbit) 
using  mid-range  weighting  parameters  (/>,,  =  170  and  =  550)  and  minimum 
control  acceleration  weighting  parameters  (/»,,  =  36  and  =  1000).  The  noise 


statistics  used  in  the  simulations  are  given  in  Table  3.1.  These  statistics  are 
consistent  with  the  tracking  accuracy  study  contained  in  the  flight  dynamics 
study  [24]. 


Table  3.1:  Simulation  Noise  Statistics 


Fixed  Power 

Average 

Spectrum 

State  a 

Position 

TTp^  =  2.665  ^ 

9.637  km 

Velocity 

7r„2  =  2.665e-8  ^ 

1.152  f 

Acceleration 

=  1.666e-15  ^ 

N/A 

The  variations  in  the  average  total  steady  state  propulsion  requirement  per  day 
(AVtot)  with  respect  to  these  input  parameters  is  summarized  in  the  next 
section  and  Appendix  B.  The  reference  input  was  a  clockwise  circular  halo 
orbit  with  a  radius  of  3500  km  and  a  non-dimensional  frequency  of  1.862647. 

Appendix  Pi.l  summarizes  the  state  .space  realization  and  the  transfer 
function  matrix  for  the  sixth-order  Hi  guidance  law  using  mid-range  parain 
eters.  The  guidance  law  is  tniniiuum  phase  and  stable.  Since  the  linearized 
equations  of  motion  are  uncoupled  between  the  ^  and  //-a.xes  and  the  (.'-axis, 
the  guidance  law  is  also  uncoupled  in  this  manner. 

Appendix  B.2  shows  the  last  10  revolutions  of  the  system  states, 
estimation  error  and  control  acceleration  when  the  system  has  reached  steady 
state.  Both  the  system  states  and  the  control  accelerations  are  sinusoidal 
functions  at  the  selected  halo  orbit  frequency.  The  estimation  error  is  quickly 
driven  to  zero  by  the  observer  in  all  states.  Table  3.2  compares  these  typical 
guidance  laws  with  previous  results;  these  results  do  not  include  proct'ss  or 
measurement  noi.se. 

.'\  rol)u.stiie.ss  analysis  was  done  for  tin*  niid-range  weight  ing  parame 


Table  3.2:  AVtot  Comparison 


Study 

AVtot  (m/s/day) 

Mid-range  Weighting 

0.2997 

Minimum  Weighting 

0.1549 

Farquhar  [13] 

0.2839 

Flight  Dynamics  [24] 

0.2800 

Heppenheimer  [18] 

0.2929 

Fraietta  and  Bond  [17] 

0.7934 

ters  guidance  law.  Structured  uncertainty  was  used  to  represent  error  .sources 
such  as  neglected  high  frequency  dynamics,  input  actuator  errors  and  low  fre¬ 
quency  plant  parameter  errors.  The  gain  and  phase  margins  were  reduced  with 
the  addition  of  the  observer  into  the  guidance  law.  The  gain  and  phase  margin 
could  be  improved  through  the  loop  transfer  recovery  technique;  however,  this 
procedure  was  not  applied  in  this  research. 

Figure  3.4  gives  the  geometric  view  of  the  halo  orbit  for  the  mid¬ 
range  weighting  parameters  'H2  optimal  guidance  law  when  both  process  and 
measurement  noise  are  absent  from  the  simulation;  Fig.  3.5  applies  when  3(7 
process  and  3(7  measurement  noise  are  included.  The  orbits  are  stable, 
repeatable  and  the  line-of-sight  with  the  Earth  is  never  lost.  Measurement 
noise  showed  little  effect  on  the  geometry  of  the  resultant  halo  orbit;  proc(*ss 
noise  showed  a  greater  impact.  Appendix  11.3  gives  additional  geumetri(  plots 
of  the  front  and  side  views  of  the  halo  orbit  with  various  amounts  of  process 
and  measurement  noise  added  to  the  simulation. 


Figure  3.1:  Halo  Orbit  (leometry  (Z<*ro  .Noise) 


Figure  3.5:  Halo  Orbit  Geometry  (3(7  Noise) 
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3.4.3  Control  Input  Parametric  Study 

Three-axes  versus  two-axes  control  and  (-axes)  was  investigated 
to  determined  the  advantages  or  disadvantages  of  each  method.  The  figures- 
of-merit  for  this  parametric  study  were  chosen  to  be  the  average  steady  state 
AV'  per  day  for  maintaining  the  orbit  ((  and  //-axes)  and  the  percent  steady 
state  error  from  the  desired  circular  halo  orbit. 

Table  3.3  gives  the  results  for  three  cases  (a-c)  where  the  semi-minor 
axis  steady  state  error  of  the  resultant  halo  orbit  was  held  constant;  the  po¬ 
sition  weighting  factors  were  varied  in  order  to  maintain  constant  semi-minor 
axis  errors.  Semi- minor  axis  errors  more  negative  than  -11.43%  indicate  a  loss 
of  line-of-sight  with  the  Earth.  In  addition,  the  average  steady  state  AV'  per 
day  in  the  (-axis  was  held  constant.  Three-axis  control  provided  a  better  re¬ 
sultant  halo  orbit  which  was  closer  to  the  desired  circular  orbit.  However,  the 
AV'  associated  with  the  tighter  control  was  two  orders  of  magnitude  higher. 


Table  3.3:  Three  vs.  Two  Axes  Control 


Control 

Axes 

<1 

Av; 

Major 

Minor 

(ni/s/day) 

(%  c 

■rror) 

a)  3 

7.866 

2.691 

4.85 

-11.. 55 

2 

0.034 

1 1 .92 

-11.  .34 

b)  3 

8.124 

2.778 

4.20 

-5.42 

2 

0.063 

8.53 

-5.26 

c)  3 

8.476 

2.890 

2.62 

-0.94 

2 

0.209 

- 

4.03 

-0.95 

Two-axes  control,  using  either  control  in  the  (  or  //-axis,  was  also 
investigated  for  advantages  or  disadvantages.  Table  3.4  gives  the  results  for 
three  cases  (a-c)  where  the  average  steady  state  AV  per  day  for  each  case  was 
held  constant;  the  position  weighting  factors  were  varied  in  order  to  maintain 
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Table  3.4:  (  vs.  ij  Axis  Control 


Control 

Axes 

S' 

<1 

Major 

Minor 

(m/s/day) 

(%e 

rror) 

a)  (f 

0.034 

- 

11.92 

-11.34 

V 

0.034 

65.84 

-27.62 

b)  e 

0.063 

8.53 

-5.26 

- 

0.063 

47.29 

-78.86 

0.209 

4.03 

-0.95 

V 

- 

0.209 

42.12 

-91.56 

r„  e„  I— —I  u 
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Figure  3.6:  Representative  Example  Block  Diagram 

constant  average  steady  state  At  per  day.  'I'wo-axes  (■ontrcd  using  f-axis 
provided  much  tighter  resultant  halo  orbits  for  the  same  propulsion  cost. 

These  two  findings  are  a  result  of  computing  the  guidance  law  using 
regulator  theory  and  then  commanding  an  input  to  be  tracked  by  the  closed- 
loop  system.  This  procedure  effectively  clo.ses  an  inner  loop  containing  the 
states  and  associated  gains  which  are  not  being  commanded  to  track;  an  outer 
loop  is  then  closed  which  contains  the  states  which  are  tracked.  Figure  3.6 
shows  these  loops  for  a  representative  linear  quadratric  regulator  example 
using  the  dynamics  of  the  halo  orbit  guidance  problem.  Figure  3.7  plots  the 
example’s  control  acceleration-to- reference  input  transfer  function  (^)  for  the 


Figure  3.7:  Control  Acceleration  Transfer  Function  (^)  for  the  Three  and 
Two  Axes  Control  Example 

three-axes  and  two-axes  control  cases.  At  the  frequency  of  interest  (1.862647), 
three-axes  control  takes  significantly  more  control  acceleration  than  either  of 
the  two-axes  control  cases.  Tracking  accuracy  is  given  by  the  magnitude  of  the 
output-to-reference  input  transfer  function  (^)  plotted  in  Fig.  3.8.  “Cood" 
command  following  for  this  example  would  uteui  al  a  magnitude  ot  0  dll. 
Hence,  three-axes  control  has  (li«*  l)est  tracking  performance  followed  by  two 
axes  ^  control;  two-axes  //  control  has  tin-  worst  tracking  performance.  .Note 
that  the  trends  shown  in  this  represeiitativr*  example  are  the  same  as  the 
trends  shown  in  Tables  3.3  and  3.4.  'I’lius,  two-axes  ^  control  is  used  in  the 


Non-dimensional  Frequency 


Figure  3.8:  Closed-Loop  Transfer  Function  {^)  for  the  Three  and  Two  Axes 

'  V 

Control  Example 


remaining  simulations  because  at  the  frequency  of  interest  it  requires  the  least 
control  acceleration  for  the  best  command  following. 

3.4.4  Noise  Characteristics  Parametric  Study 

The  steady  state  solution  to  the  observer  Ricatti  equation  i.s  the 
error  covariance  associated  with  the  estimated  system  states.  Equation  (3.19) 
shows  the  error  covariance  is  a  function  of  77^  and  The  noise  statistics 
used  in  the  simulations  are  given  in  Table  3.1.  These  statistics  arc  consistent 
with  the  tracking  accuracy  study  contained  in  the  flight  dynamics  study  [24]. 
Figu'^e  B.6  (a)  in  Appendix  B.4  shows  the  average  standard  deviation  of  the 
position  and  velocity  state  as  a  function  of  the  standard  deviation  of  the 
process  noise;  Fig.  B.6  (b)  gives  the  average  standard  deviation  of  the  velocity 
state.  In  both  figures,  the  fixed  power  spectrum  of  the  measurement  noise  is 
assumed  to  be  2.665  m^/s  in  each  axis.  The  average  standard  deviation  of 
the  position  state  is  defined  as  the  numerical  average  of  th('  position  singular 
values  of  the  error  covariance  matrix.  The  similar  definition  appli('s  to  the 
average  standard  deviation  of  the  velocity  state. 

A  monte  carlo  analysis  was  completed  to  quantify  the  effect  of  these 
noise  statistics  on  the  average  steady  state  AF  per  day  for  halo  orbit  mainte¬ 
nance.  Thirty  150  day  simulations  were  run  and  the  resultant  average  steady 
state  AF  per  day  obtained.  In  all  cases,  the  reference  input  to  the  simula¬ 
tion  was  a  clockwise  circular  halo  orbit  with  a  radius  of  3500  km  and  a  non- 
dimensional  frequency  of  1.862647.  Mid-range  weighting  parameters  were  also 
used.  Figure  3.9  shows  essentially  a  linear  relationship  between  the  amount  of 
process  and  measurement  noise  and  the  resultant  average  steady  state  AV’  per 
day.  The  dashed  lines  show  the  minimum  and  maximum  values  obtained  in 
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the  30  simulations.  The  average  steady  state  Al/  per  day  found  in  Table  3.2 
corresponds  to  zero  process  noise  and  zero  measurement  noise. 

3.4.5  Halo  Orbit  Characteristics  Parametric  Study 

The  position  weighting  factors,  p,,  and  and  the  halo  orbit  fre¬ 
quency  were  varied  to  determined  the  effects  of  each.  Non-dimensional  fre¬ 
quency  was  varied  from  0.2  to  7.4  for  each  case.  As  before,  the  reference  input 
to  the  simulation  was  a  clockwise  circular  halo  orbit  with  a  radius  of  3500  km. 
The  duration  of  each  simulation  was  150  days. 

Appendix  B.5  gives  the  results  of  this  parametric  study  as  a  func¬ 
tion  of  halo  orbit  frequency.  In  each  case,  the  resultant  steady  state  halo 
orbit  was  an  ellipse.  For  halo  orbits  with  clockwise  orbital  rotation.  Fig.  3.10 
shows,  qualitatively,  the  effect  on  the  halo  orbit  of  increasing  the  desired 
halo  orbit  frequency.  The  minimum  deviation  from  a  purely  vertical  orien¬ 
tation  was  obtained  at  a  non-dimensional  frequency  of  1.17.  As  the  halo  or¬ 
bit  non-dimensional  frequency  deviated  from  the  system's  iialiiral  frccjueiic  y 
(1.862647),  velocity  requirements  in  both  a.xes  increa-sed  rapidl}.  When  non 
dimensional  frequency  was  reduced  by  1,  the  average  steady  state  AV'  per  day 
increased  by  a  factor  of  30.  For  an  increase  in  non-dimensional  frequency  of 
1,  the  average  steady  state  AV  per  day  increased  by  a  factor  of  50. 

The  halo  orbit  non-dimensional  frequency  was  then  fixed  at  the  nat¬ 
ural  frequency  of  the  system  and  the  position  weighting  factors  were  varied; 
these  results  are  shown  in  Appendix  B.6.  The  geometry  of  the  halo  orbit 
was  essentially  constant  with  respect  to  changes  in  the  (-axis  weighting  fac¬ 
tor  and  steady  state  error  decreased  as  the  ?/-axis  weighting  factor  increased. 
The  orientation  of  the  halo  orbit  was  essentially  constant  for  changes  in  either 


\  (Non-dimensional) 
(b)  Side  View 


Figure  3.10:  Effect  of  Increasing  Frequency  on  Halo  Orbits  with  Clockwise 
Orbital  Rotation 


weighting  factor.  Due  to  the  uncoupling  of  the  linearized  equations  of  motion, 
the  velocity  required  in  the  ^-axis  Wcis  only  influenced  by  the  weighting  factor 
in  the  T/-axis;  the  velocity  required  increased  as  the  weighting  factor  increased. 
Similarly,  the  velocity  required  in  the  ^-axis  was  only  influenced  by  the  weight¬ 
ing  factor  in  the  (-axis.  However,  in  this  case,  the  velocity  required  decreased 
as  the  weighting  factor  increased.  This  was  due  to  the  (-axis  originally  being 
stable  and  the  (  —  »/  axes  having  to  be  stabilized  by  the  guidance  law. 


Chapter  4 


Rendezvous  Targeting 


Terminal  phase  rendezvous  for  the  Ceiniiii.  Apollo  and  Space  Shutt  le 
missions  was  executed  from  stable  elliptical  orbits  in  the  two- body  problem. 
Two-body  motion  is  planar;  the  rendezvous  profile  contains  zero,  or  at  least 
small  amplitude  sinusoidal,  out-of-plane  relative  motion.  Rendezvous  target¬ 
ing  was  done  using  Lambert’s  Theorem.  By  contrast,  motion  in  the  vicinity 
of  the  translunar  equilibrium  point  is  three-dimensional  and  unbounded.  Fur¬ 
ther,  Lambert’s  Theorem  does  not  apply  to  motion  in  the  three-body  problem. 

This  chapter  considers  the  terminal  phase  rendezvous  of  two  space¬ 
craft  about  the  translunar  equilibrium  point.  It  assumes  the  target  vehicle  is 
in  a  small  radius  halo  orbit  and  is  controlling  its  own  orbit  using  an  unknown 
guidance  law;  the  target  vehicle  remains  passive  with  respect  to  the  rendezvous 
profile.  The  chaser  vehicle  is  in  an  independent  small  radius  halo  i>rldt  and 
executes  all  rendezvous  maneiiv<‘rs.  I'liis  <hapt<‘r  l)egins  In  (le\(‘loping  tin* 
state  equations  for  the  target  and  chaser  vehicles,  rin*  targeting  law  is  then 
derived;  this  derivation  is  loosely  equivalent  to  the  rendezvous  application  of 
Hill’s  equations  in  the  two-body  problem.  Lastly,  nonlinear  simulation  results 
are  presented. 

4.1  State  Equations 

The  target  vehicle  orbits  de.srribed  in  ('liapter  3  can  be  approximated 
by  an  ellipse  inclined  to  the  r/  -  (  plane.  The  equations  of  motion  for  this 
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approximation  are 


'^iV  = 


^(r7iax  “i"  ) 

Vtmax  sin(ft  +  I/.',,) 
C(maxSin(ft  +  ll'^) 


and 


which  leads  to  the  state  equation 


K 

—  A 

'^tp 

-  Ap 

(i.l) 


(4.2) 


where 


(4.3) 


and  f  is  the  halo  orbit  frequency.  The  subscripts  (  ),p  and  (  )<„  denote  the 
target  vehicle’s  position  and  velocity,  respectively. 


Since  the  halo  orbits  being  considered  have  a  small  radius,  the  lin¬ 
earized  equations  of  motion,  Kq.  (2.40).  for  the  cha.ser  vehicle'  can  be  used. 
Further,  during  a  rendezvous  profile  the  chaser  vehicle  will  t('riiiiuate  its  own 
halo  orbit  guidance.  Hence,  the  chaser  vehicle's  state*  e'eiuatie)n  be'eoiue's 


^cp 

^cv 


=  Aa 


/?C 


cp 


(4.4) 


where  Ac  is  defined  by  Eq.  (2.43)  and  the  subscripts  (  )fp  and  (  )ct,  denote 
the  chaser  vehicle’s  position  and  velocity,  respectively. 

Defining  the  relative  position  and  velocity  as 


'6p 

A 

Ocp  -  ^,p 

Oc.  - 

leads  to  the  relative  state  equation 


6p  ] 

A(; 

( Ac;  -  A,, ) 

[  ^  1 

6, 

*tp 

.  . 

0 

A, 
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Equation  (4.6)  can  be  solved  analytically 


«p(^2-)  ^ 

^11  ^13  ‘Ah 

r 

^v(^2  ) 

<f>2i  4>22  ^23  ^24 

«v(<l  +  ) 

0  0  <^33  <^34 

iPtv(h-)  . 

.  0  0  043  044  - 

where  ii  and  I2  denote  the  initiation  and  final  limes  of  the  rendezvous,  re¬ 
spectively.  Also,  the  state  transition  matrix  is  a  function  of  ('xponentially 
increasing  and  decreasing  terms  as  well  as  sinusoidal  terms  at  three  different 
natural  frequencies;  is  further  detailed  in  Appendix  A. 2.  The  super¬ 
scripts  (  )~  and  (  )■*■  denote  immediately  before  and  immediately  after  a 

specified  time,  respectively. 


4.2  Targeting  Law 

The  necessary  and  sufficient  conditions  for  rendezvous  are  that  the 
relative  position  and  velocity  at  the  time  of  rendezvous  be  a  specified  value. 
Hence,  and  ^,,(<2'*’)  are  known  paramet<*rs.  (’onsider  a  two  maneuver 

terminal  pha.se  rendezvous  seqiKuice  where  the  first  maneuver  is  exec  uled  at 
the  beginning  of  the  rendezvous  (f,)  and  designed  so  that  the  chascu  vcdiic  le 
will  coast  to  the  specified  relative  position  at  the  time  of  reudt'zvous  (/j); 
the  second  maneuver  is  then  executed  to  satisfy  the  s|)ecilied  relative  vcdocily 
condition. 

The  terminal  pha.se  initiation  maneuver  (T/)  is  defined  a.s 

avt-/  =  -  aMr)  =  -  «.(<.')•  (4.8) 

Given  the  first  row  of  the  state  transition  matrix  from  Kq.  (4.7),  Eq.  (4.8)  can 
be  expressed  in  terms  of  the  specifiecl  position  at  the  end  of  the  rendezvous;  the 


initial  relative  position  and  velocity;  and  the  chaser  vehicle’s  initial  position 
and  velocity 

AVtI  =  4>i2  '  [^p(^2^)  +  {013  “  )  +  {0H  -  ) 

-0i3ap(<r)-0H^>c.(<r)).  (1.9) 

Note  that  this  targeting  law  requires  only  relative  state  information,  which 
can  be  provided  by  an  onboard  rendezvous  navigation  filter  using  proximity 
sensors,  and  chaser  vehicle  state  information.  Hence,  no  target  vehicle  spe¬ 
cific  information  is  required.  Additionally,  Eq.  (4.9)  contains  a  numerical 
singularity  at  specific  transfer  times.  That  is,  when 

(t2  -u)=  -  k  7.637724  days;  1,2,3, •••  (4.10) 

the  <^129  element  becomes  zero  causing  the  0,.2  submatrix  to  become  singu¬ 
lar.  Thus,  0|2~'  does  not  exist.  The  <i>i29  element  corr(*sponds  to  the  un¬ 
coupled  (-axis  in  the  lineariz<‘d  eciuations  of  motion.  I'lie  ( compoiKMit  is  a 
sinusoidal  function  and  Ix'comes  zero  twice  during  a  cy<le  at  its  natural  fre 
quency  (1.786178).  These  singularitit's  do  not  limit  this  |•(*s(‘arch  bi'cau.se  the 
transfer  times  considered  here  are  much  sliort<*r  than  any  singularity  time. 

Lastly,  the  terminal  phase  final  maneuver  (I  F  )  is  defined  as 

AVrr  =  «cv(<2-')  -  «c.(<2-)  =  6M2^)  -  (4.11) 


4.3  Simulation  Results 

A  nonlinear  simulatioti  wa.s  developed  to  model  the  three-body  ter¬ 
minal  phase  rendezvous  problem  (stn*  Chapter  7).  for  analysis  purposes,  an 
option  exists  which  will  place  the  target  vehicle  in  a  planar  circular  orbit  in 
dined  to  the  r/-(  plane  by  a  specified  incline  angh*  rather  than  using  an  actual 
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To  Earth 


T1 


Figure  4.1:  Incline  and  Initial  Condition  Angle  Definitions 

halo  orbit.  This  option  eliminates  the  target  vehicle  trajectory  approximation 
used  in  the  development  of  Eq.  (1.9).  In  addition,  an  initial  condition  angle 
was  defined  to  indicate  the  position  on  the  halo  orbit  of  the  T I  maneuver. 
Figure  4.1  shows  the  geometry  of  these  angles.  Several  profiles  were  executed 
in  order  to  investigate  the  characteristics  of  the  three- body  rendezvous  prob¬ 
lem.  For  each  profile,  clockwise  orbits  were  used.  The  target(*(l  position  and 
velocity  offset  at  rendezvous  was  specified  to  be  z«to  in  all  t  as«  s. 

4.3.1  TI  Targeting  Law  Error  Parametric  Study 

The  first  three-body  rendezvous  characteristic  investigated  was  final 
position  error  at  the  time  of  rendezvous  due  to  the  77  targeting  law.  The 
targeting  la  in  Eq.  (4.9)  leads  to  a  certain  amount  of  error  at  the  time  of 
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rendezvous  due  to  the  target  vehicle’s  trajectory  approximation  introduced 
in  its  development.  Errors  from  the  linearization  of  the  equations  of  motion 
for  the  chaser  vehicle  have  been  removed  by  iterating  on  the  specified  relative 
position  at  the  time  of  rendezvous  within  the  simulation;  this  type  of  iteration 
technique  is  also  used  in  the  Space  Shuttle  rendezvous  targeting  software. 

Figure  4.2  (a)  gives  the  results  for  a  clockwise  halo  orbit  as  a  function 
of  transfer  time  and  initial  condition  angle.  Figure  1.2  (b)  show.s  couslant 
error  contours  from  f'ig.  4.2  (a).  Figure  4.2  (c)  plots  the  side  view  of  the 
target  vehicle  trajectory  using  the  inclined  ellipse  approximation  (dotted  line) 
with  the  side  view  of  the  actual  halo  orbit  (solid  line).  The  maximum  final 
position  error  occurs  when  T I  occurs  at  the  left  or  right  sides  of  the  halo 
orbit  (initial  condition  angle  approximately  90  or  270  deg).  Further,  the  final 
position  error  was  generally  in  the  out-of-plane  direction.  Figure  4.2  (c)  show 
that  these  points  of  maximum  error  also  corresponds  to  the  points  of  maximum 
curvature  of  the  actual  halo  orbit.  For  all  profiles,  the  chaser  vehicle  initiated 
the  rendezvous  profile  approximately  19  km  behind  the  target  vehicle. 

To  correct  for  this  targeting  law  error,  a  tnidcoursc'  correction  ma¬ 
neuver  was  added  to  the  rendezvous  profile.  For  deinoiist ration  purposes,  the 
transfer  lime  was  held  constant  at  9  hrs  and  llu'  initial  (•i)nditii)n  aiigh-  was 
fixed  at  46.08  deg;  the  transfer  time  and  initial  comlilion  angle  were  arl)ilrarily 
chosen  to  be  these  values  to  be  consistent  with  the  planar  circular  orbit  ba.se- 
line  case  that  results  in  a  completely  “in-plane”  rendezvous  profile.  Figure  4.8 
shows  the  results  for  one  midcourse  correction  maneuver  executed  between 
TI  and  TF  minus  one  hour.  I'lie  final  position  error  decreast'd  rapidly  from 
1122  km  to  18  km  while  the  total  propulsion  cost  of  the  rendezvous  increa.sed 
only  slightly  from  1.25  ni/s  to  1.78  ni/s.  Hence,  the  addition  of  niidcourse  cor- 
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Figure  4.3:  Midcourse  Correction  Maneuver  Placement 

rection  maneuvers  can  substantially  reduce  errors  due  to  the  target  vehicle’s 
trajectory  appro.ximation  used  in  the  development  of  Eq.  (4.9).  The  optimal 
location  of  a  midcourse  correction  maneuver  can  not  be  determined  solely 
from  these  data  since  that  analysis  requires  a  rendezvous  navigation  filter  to 
be  present. 

4.3.2  Constant  Incline  Angle  Parametric  Study 

The  .second  three-body  rendezvous  characteristic  investigated  was 
the  effect  of  transfer  time  and  initial  condition  angle  on  the  geometry  of  the 
rendezvous  profile.  A  planar  circular  orbit  was  used  to  insure  final  position 
errors  were  less  than  one  meter.  The  incline  angle  was  fixed  at  22.85  deg  whicli 
corresponded  to  the  incline  angle  of  the  actual  halo  orbit  u.sed  in  this  research. 
Finally,  the  chaser  vehicle  initiated  the  rendezvous  profile  19  km  behind  the 
target  vehicle  in  a  coelliptic  orbit. 

Figure  4.4  (a)  shows  the  niaxiiniini  out-of-plane  relative  position 


achieved  by  the  rendezvous  profile  as  a  function  of  the  input  parameters. 
The  out-of-plane  component  results  from  the  coupling  of  the  ^  —  t)  axes  in  the 
equations  of  motion.  Its  magnitude  is  greatest  when  the  rendezvous  profile 
is  initiated  at  the  top  or  bottom  of  the  orbit  (initial  condition  angle  approx¬ 
imately  0  or  180  deg).  These  points  correspond  to  the  points  of  maximum 
distance  from  the  equilibrium  point  along  the  action  line  of  gravitational  and 
centrifugal  forces.  Figure  4.4  (b)  gives  the  total  propulsion  cost  of  the  ren¬ 
dezvous  profile;  Fig.  4.4  (c)  overlays  a  plot  of  a  minimum  out-of-plane  ren¬ 
dezvous  profile  (dashed  lines)  on  constant  total  propulsion  cost  contours  (solid 
lines).  Two  characteristics  can  be  observed  from  these  plots.  First,  for  every 
initial  condition  aiigle,  there  is  a  transfer  time  which  produce  a  miiiiiinnii  total 
propulsion  cost  of  the  rendezvous  profile  (eg.  180  deg,  18  hrs).  .Second,  total 
propulsion  costs  increase  with  out-of-plane  relative  motion. 

4.3.3  Constant  Transfer  Time  Parametric  Study 

The  third  three- body  rendezvous  characteristic  investigated  was  the 
effect  of  incline  angle  and  initial  condition  angle  on  the  geometry  of  the  ren¬ 
dezvous  profile.  As  before,  a  planar  circular  orbit  was  used  to  insure  final 
position  errors  were  less  than  one  meter.  The  transfer  time  was  fixed  at  the 
baseline  ca^e  value  of  9  hrs  and,  once  again,  the  cha.ser  vehicle  initiated  the 
rendezvous  profile  19  km  behind  the  target  vehicle  in  a  coelliptic  orbit. 

Figure  4.5  shows  the  same  general  information  as  Fig.  4.1  except 
incline  angle  has  replaced  transfer  tim**.  I'liesc*  r**sults  show  that  a  mini 
mum  total  propulsion  cost  rendezvous  profile  is  a  function  of  incline  angle  as 
well  as  initial  condition  angle  and  transfer  time.  Also,  when  transfer  time 
is  held  constant,  the  minimum  total  propulsion  cost  rendezvous  profiles  have 
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“in-plane”  relative  motion. 

4.3.4  Representative  Relative  Trajectories 

Figure  4.6  summaries  the  geometry  of  ihe  rendezvous  profile  as  ini¬ 
tial  condition  and  incline  angles  change.  Figure  4.6  (c)  gives  a  completely 
“in-plane”  profile  executed  on  the  left  side  of  the  orbit.  Figures  4.6  (a)  and 
(b)  show  that  for  a  fixed  incline  angle  and  increasing  initial  condition  angle 
in  the  vicinity  of  the  “in-plane”  profile,  the  out-of-plane  component  switches 
from  away  from  the  Moon  (positive  to  toward  the  Moon  (negative  ^).  Fig¬ 
ures  4.6  (d)  and  (e)  show  this  same  characteristic  for  a  fixed  initial  condition 
angle  and  decreasing  incline  angle  in  the  vicinity  of  the  “in-plane”  profile. 

Figure  4.6  also  shows  that  the  relative  trajectory  is  three-dimensional. 
To  further  examine  this  conclusion,  the  torsion  and  rate  of  change  of  torsion 
was  computed  for  the  relative  trajectory  at  the  initiation  of  the  rendezvous 
profile  (Tl).  Torsion  is  a  measure  of  the  twisting  of  a  (•urv(>  and  is  (()mput(‘ 


r  =  (4.12) 

!l  6,  X  6,  IP 

In  addition,  the  rate  of  change  of  torsion  is  given  by 

j.  _  (^p  X  6^.)'  6p  d~(^p  X  Sp)'  Sp  —2t(Sp  X  6p)  •  (6p  x  6p) 

ll«pX«p!P 

Figures  4.7  (a)  and  (b)  plot  the  relative  trajectory’s  torsion  and  rate 
of  change  of  torsion  at  the  beginning  of  the  terminal  phase  as  a  function  of 
initial  condition  and  incline  angles.  Transfer  time  was  held  constant  at  9 
hrs.  A  necessary  condition  for  the  relative  trajectory  to  be  two-dimensional 
is  that  the  torsion  and  rate-of-change  of  torsion  be  zero  at  the  initiation  of 
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Figure  4.6:  Representative  Relative  Motion 
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the  rendezvous  profile.  Figure  4.7  (c)  shows  contours  of  torsion  and  the  rate 
of  change  of  torsion  equal  to  zero.  Also  shown  is  the  minimum  out-of-plaiie 
contour  from  Fig.  4.5.  Figures  4.7  (d)  and  (e)  are  expanded  views  of  the  inter¬ 
sections  of  the  three  contours  given  in  Fig.  4.7  (c).  For  the  relative  trajectory 
to  be  two-dimensional  and  lie  in  the  target  vehicle’s  halo  orbit  plane,  the 
three  contours  must  intersection  at  the  same  point.  Hence,  two-dimensional 
relative  trajectories  which  lie  in  the  target  vehicle’s  halo  orbit  plane  do  not 
exist.  However,  the  minimum  out-of-plane  contour  from  Fig.  4.5  represents 
a  maximum  out-of-plane  distance  on  the  order  of  10  m.  Thus,  thesr*  relative 
trajectories  are  essentially  “in-plane”  and  will  be  referred  to  as  such  in  this 
research. 


Chapter  5 


Rendezvous  Navigation 


The  Apollo  onboard  rendezvous  navigation  filter  represents  the  first 
application  of  Kalman  filtering  techniques  to  the  rendezvous  navigation  prob¬ 
lem.  This  problem  can  be  defined  as  determining  th(‘  relative  staU'  of  th«'  two 
spacecrafts  performing  the  rendezvous.  .Muller  and  Kachmar  [10]  showed  tlu' 
equivalence  of  optimally  estiinaliiig  the  relative  stale  or  optimally  estimating 
the  inertial  state  of  one  spacecraft  and  subtracting  it  from  an  externally  sup¬ 
plied  inertial  state  of  the  other  spacecraft.  This  estimation  is  optimal  in  the 
sense  that  it  yields  the  minimum  relative  state  error  at  the  time  of  measure¬ 
ment  incorporation,  not  necessarily  the  minimum  relative  state  error  over  the 
entire  navigation  period.  Nevertheless,  to  reduce  the  complexity  of  the  ren¬ 
dezvous  navigation  filter,  this  formulation  was  used  in  the  Apollo  and  Space 
Shuttle  rendezvous  navigation  filters. 

For  this  research,  the  rendezvous  navigation  filter  will  estimate  the 
inertial  state  of  the  chaser  vehicle;  it  is  assnme<l  t  hat  either  the  target  vehicle  or 
its  mission  control  center  will  supply  the<  ha,s<*r  vehicle  with  target  vc'hicle  stale 
information  at  the  time  of  each  measiireimuit.  I  Ik’  Space  Shuttle  rendezvous 
navigation  filter  is  mechanized  in  this  mamier  except  only  one  target  vehi(  le 
state  vector  is  required  at  the  beginning  of  the  rendezvous  profile.  Target 
vehicle  state  information  is  required  at  every  measurement  in  this  research 
because  the  target  vehicle  requires  active  guidance  to  maintain  its  own  halo 
orbit  and  it  is  assumed  that  this  guidance  law  is  unknown  to  the  chaser  vehicle. 
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The  chaser  vehicle  would  have  to  know  the  target  vehicle’s  guidance  law  in 
order  to  propagate  target  vehicle  state  information  in  a  manner  similar  to  the 
Space  Shuttle  rendezvous  navigation  filter. 

In  this  chapter,  a  continuous  dynamics,  discrete  measurement  ex¬ 
tended  Kalman  filter  is  developed.  The  rendezvous  navigation  filter  is  for¬ 
mulated  assuming  the  estimated  states  are  the  chaser  vehicle's  inertial  stale, 
unmodeled  accelerations  on  the  chaser  vehicle,  and  measurement  biases.  I'he 
filter  is  implemented  in  the  nonlinear  sinuilatioii  such  that  the  cha.sei  vehicle's 
inertial  state  is  always  estimated  and  the  other  estimation  states  are  selec  table. 
Initially,  the  dynamical  system  is  developed.  Then,  the  estimation  state-error 
covariance  matrix  propagation  and  update  cKjuations  are  presented  and  the 
measurement-state  relationships  derived.  Lastly,  nonlinear  simulation  results 
are  given  to  assess  the  performance  of  the  rendezvous  navigation  filter. 

5.1  Continuous-Discrete  Kalman  Filter 

5.1.1  Dynamical  System 

The  dynamical  system  for  this  problem  can  be*  described  by  a  con¬ 
tinuous  nonlinear  differential  equation 

X,{t]  =  fi(XAt))  +  a(l)  (0.1) 

where  /,(Xc(0)  ’s  given  by  Eq.  (2.27)  and  a(t)  rc'prc'sents  immodeled  accel¬ 
erations.  Assume  the  unmodeled  accelerations  have  the  form  of  a  first-order 
Gauss- Markov  process 

d{t)  =  Aaa(t)  ^  Wa{i)  (5.2) 


where 


j4a  —  diag(aall,aa22,fla33)aa-H.«ii55.au66] 


(5.3) 


64 


and  Wa{t)  is  white  gaussian  noise  with  zero  mean  and  fixed  power  spec¬ 
trum.  The  discrete  nonlinear  measurements  are  modeled  a.s 

y{tk)  =  /i,(X,(4))  +  b{lk)  +  v(4);  ^  =  1,2,  •  •  •  (5.4) 

where  represents  the  measurement- state  relationships  derived  in 

Section  5.1.3.  The  mezisurement  bias  6(4)  is  also  modeled  as  a  first-order 
Gauss- Markov  process 

b{t)  ^  Abb{t) Wh{t)  (5.5) 

where 

Ab  =  diag  [«(,!!,  <444  (5.6) 

and  W(,(l)  is  white  gaussian  noise  with  zero  mean  and  Q/^  fixed  power  spectrum. 
The  measurement  noise  t>(4)  is  also  white  gaussian  noise  with  zero  mean  and 
Ri  covariance  matrix. 

Define  an  estimation  state  vector 

Z(l)=  0(0  .  (5.7) 

6(0 

Combining  Eqs.  (5.1),  (5.2),  and  (5.5)  yields 

Z(t)  =  (5.8) 

which  for  notational  conveni<Tice  will  be  written  as 

Z  =  f  2-^  w  (5.!)) 

where 

r/,(X.(<))  +  a(<)' 

/2  =  AMi)  (5.10) 

AMt) 

0  ■ 

w  =  ti;(/)  =  Wait)  . 

Wh{t)  , 


(5.11) 


65 


The  white  gaussian  noise  w  has  zero  mean  and  Q  fixed  power  spectrum  where 


Q  = 


0  0  0 

0  Qa  0 

[0  0  Q,\ 


(5.12) 


Finally,  Eq.  (5.4)  becomes 


Vih]  =  hj(.Z((t))  +  »(U) 


(5.13) 


which  for  notational  convenience  will  be  written  as 


Vk  =  +  Vk 


(5.14) 


where 


h2k  =  h,{Z{ik))  =  fi,(A',(4))  +  6(4)  =  hyk  +  6,.  (5.15) 


5.1.2  Propagation  and  Update  Equations 

Equations  (5.9)  and  (5.14)  give  the  standard  form  of  the  dynamical 
system  for  the  continuous-discrete  Kalman  filter.  The  derivation  of  the  prop¬ 
agation  and  update  equations  for  this  filter  is  well-known  (see  Gelb  [41]  for 
example);  this  derivation  minimizes  the  magnitude  of  the  expected  value  of 
the  estimation  error.  The  estimation  state  propagation  equation  is 

^  =  E(/,l.  (5.16) 

Expanding  about  the  estimated  state  in  a  Taylor  series,  taking  the  expec¬ 
tation  and  neglecting  higher  order  terms  yields 


z  =  f,(zm 


4  <  /  <  4-n 
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where  sensor  measurements  are  available  at  If,  and  lk+\.  For  iiolat ioiial  ((Ui- 
venience,  Eq.  (5.17)  will  be  written  as 

i  =  /j.  (5.18) 


Propagation  of  the  error  covariance  matrix  P{t)  is  done  using 

P  =  PF^  +  FP  +  Q 


(5.19) 


where  Q  has  the  form  given  by  Eq.  (5.12)  and 


dfmt)) 


dZ(l) 


Equation  (5.20)  can  be  further  reduced  to 

■  I  ^  ° 

dX.{t)  \xAt)=X.it) 

0  0 
0  0  Ab  . 

which  after  compacting  the  notation  becomes 


L  0  0  Ab  J 

Differentiating  /,,  given  in  Fiq.  (2.27),  yields 


■  0 

0 

0 

i 

0 

0 

0 

0 

0 

0 

1 

0 

_ 

0 

0 

0 

0 

0 

1 

dx. 

Ui 

U 

0 

2 

0 

f.u 

fs2 

hs 

_2 

0 

0 

h. 

0 

0 

0 

(5.20) 


(5.21) 


(5.22) 


(5.2;n 
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where 


/41 


/42 

/43 

/52 

/sa 

/ea 

PMC 

Pec 


j  ^  3p(x,  -  (1  -/i))^  _  P 

PMC^  PA/C^ 

I  3(1  -p)(jc  +  p)^  _  (1  -  P) 

Pec^  Pec^ 

3p(xc  -  (1  -  p))Pc  ^  3(1  -  /0(.re  +  n)yc 
Pmc^  Pec’’ 

3p(Jc  -  (1  -  p))^.  3(1  -  /0(.r,  -t-  p)c 

PA/r^  Pec’ 

1  +  ^  ^  3(1  -  _  (1  -  //) 

Pa/c®  Pmc^  Pec^  Pec'^ 

^pilcZc  ^  3(1  -  p)pcic 
Pa/c®  Pec^ 

3pic^  _  P  3(1  -p)ic^  _  (1  -  p) 
PAfC^  PAfC^  Pec^  Pec^ 

\/(ic  -  (1  -  p))^  + 

+  p)2  +  +  Zc^. 


(5.24) 

(5.25) 
(  ■')  26) 


(5.28) 

(5.29) 

(5.30) 

(5.31) 


The  (  )c  subscript  denotes  the  chaser  vehicle.  Finally,  the  estimation  state 
and  error  covariance  matrix  are  updated  according  to 

Z+  =  Z-  +K,(y,-hn)  (5..32) 

=  (I  -  K,H,)P-(I  -  K,H,y  +  K,R,K,‘  (5.33) 


where  the  superscripts  (  )”  and  (  )■*■  denote  iuiiiiediately  before  and  iriinie- 
diately  after  measurement  incorporation,  respectively,  and 


K, 

(5.34) 

^2* 

=  E(ft2*] 

(5.35) 

H, 

A  dh^iZit,)) 

Z(/t)=Z(<,-)’ 

(5.36) 
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Expanding  /12  about  the  estimated  state  in  a  Taylor  series,  taking  the  expec¬ 
tation  and  neglecting  the  higher  order  terms  in  Eq.  (5.35)  yields 

hn  =  h2(Z{t,-)).  (5.37) 


Additionally,  Eq.  (5.36)  can  be  further  reduced  to 
Hk  = 


dh,(XM) 

0  I 

RxMk) 

XAhhX.(ir) 

which  after  compacting  the  notation  becornt's 

r  dh 


H,  = 


dXr 


0  I 


■ 

In  the  next  section,  h\k  is  defined  and  is  computed. 

C 


(5.38) 


5.39) 


5  1.3  Measurement-State  Relationships 

Consider  a  set  of  measurements  consisting  of 

■  R  ■ 

h  ~ 

.  0  . 


(5.10) 


where  R  is  the  relative  range  between  tlu*  chaser  and  laig«‘l  \«-lii(  les,  R  is  the 
relative  range  rate,  and  0  and  0  are  th<*  roll  and  pitch  angles  of  I  lie  relative 
position  vector  in  the  chaser  vehicle’s  body  coordinates.  I'lie  corresponding 
measurement  bias  vector  is 

{  bn  ^ 


bk  = 


(5.41) 
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Differentiating  Eq.  (5.40)  yields 


dhxk 

dXc 


r  dR 


dXc 

dk 


L  dXc  J 


(5.42) 


Assume  the  proximity  sensor  is  located  at  a  known  position  (6s)  in 
the  chaser  vehicle’s  body  coordinate  system.  Further,  assume  the  chaser  ve¬ 
hicle’s  attitude  with  respect  to  the  local  vertical  inertial  rectangular  (LVIR) 
coordinate  system  and  the  orientation  of  the  sensor's  roordinat('  system  are 
known.  Hence,  the  transformations  from  LVIH  to  the  body  coordinate  systen; 
(Tbl)  and  from  the  body  to  the  sensor  coordinate  system  (Tsh)  can  b('  com¬ 
puted.  Finally,  with  the  chaser  vehicle’s  inertial  position  and  vcdocity  known, 
the  transformation  from  the  inertial  to  the  LVIR  coordinate  system  (T/,/)  can 
also  be  computed. 


If  the  state  vector  of  the  chaser  and  target  vehicles  is  partitioned 
into  inertial  position  and  velocity, 


X.  = 

X,  = 


*•<- 

Vc 

rt 

Vi 


(5.43) 

(5.44) 


then,  the  relative  position  and  velocity  vectors  from  the  sensor  location  to  the 
target  vehicle  can  be  written 


R  =  {Thi.Tjj)'  6s 


V 


Vt  -  V^.. 


(5.45) 

(5.46) 


From  Eqs.  (5.45)  and  (5.46),  it  follows  that  the  relative  range  and  range  rate 
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Figure  5.1:  Angle  Measurements 


Further,  differentiating  Eq.  (5.46)  gives 


dV 

dX, 


(5.57) 


Figure  5.1  shows  the  two  angles  of  the  relative  position  vector  mea¬ 
sured  by  the  sensor.  These  angles  are  in  the  sensor  coordinate  frame  and  ar«‘ 
computed  by 


Q  =  arctan  [  ^  j  (5.58) 

\R-ks) 

0  =  arctan  (  .  R  is  |  (5.59) 

\sJ(Rhy+{R-ksy)l 


Figure  5.2  gives  the  two  angles  used  by  the  extended  Kalman  filter;  they  are 
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A 


Figure  5.2:  Computed  Angles 


the  preferred  line-of-sight  angles  because  they  have  an  intuitive  interpretation 
to  the  crewman  on-board  the  chaser  vehicle.  These  angles  are  in  the  chaser 
vehicle’s  body  coordinate  system  and  can  be  computed  directly  from  the  sensor 
measured  angles  by 


<i> 

8 


arcsin 


f( 


Tsb{  1 , 2)  sin  ^  2)  sin  a  cos  0 

+T’sb(3,  2)  cos  a  cos  /J 


arcsin 


i„( 


cos  8 

V  / 

-Tsb(  1 , 1 )  sin  -I-  T  ss(2,  1 )  sin  q  cos  (3 

—  TsB{^y  l)cOSQCOS/^ 


f 

] 


(5.60) 

(5.61) 
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Differentiating  Eqs.  (5.60)  and  (5.61)  yield 
d(i>  sin  ^  sin  0  d6 


dX, 


cos  ^  cos  6  dXc 

T5fl(2, 2)cosdcos/8  +  T5B(3,2)sindcos^  da 
cos  cos  0  dX, 

I  Tsb(  1,2)  cos  4  +  r5B{2, 2)  sin  d  sin /) 

— Tss(3,2)cosdsin  ft 


+ 


dft 


de 

dXc 


cos  ^  cos  ^  dXc 

^sb(2,  1)  cos  a  cos  ft  +  1)  sin  d  cos/)  da 


cos  9 

^sb(I)  1)cos4  +  Tsb{^ 

-Tsb(3,  l)cos 


2, 1 )  sin  d  sin /)  \ 
a  sin  ft  ) 


dXc 

dft 


cos  9 


Further,  differentiating  Eqs.  (5.58)  and  (5.59)  yield 


dX, 


da 


dft 

dX, 


+  (R^k 

R^ks 


R^ks)^  \  dX, 


+  ks^ 


dR  ' 

dX,, 


iRTjsy  +  {RTks)A  dXc  dxj 

(r 


+  (fl^s)^  (i,T  dh  ^  J^r  dR 


+  ■ 


(R^is)(R^  ks) 


dX, 


__  ( A’ 

tv)'  V 


dx, 


+  ib,s- 


,  dR 

dX 


:5.62) 


(5.63) 


(5.64) 


+ 


\Rip^/{¥isF+lWi 

— pjm±K.  _  f  AT^  ) ,  ,,.65 

||A||V(A''js)'  +  (A^'ts)'  V  dX,  ■  dXj 


5) 


Note  that  the  *5,  J5,  and  ks  unit  vectors  are  the  columns  of  the  sensor  to 
inertial  coordinate  system  transformation  matrix 


[*5  i  is  i  ks]=iTsBTBLTuf.  (5.66) 

Define  three  vectors 

[Cj  I  C3  i  cj  ={T.s-«Tflt)^  (.5.67) 
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which  are  not  a  function  of  Xc-  Then, 


is 

=  Tu^C-i 

(5.68) 

is 

=  Tu^C^ 

(5.69) 

ks 

=  Tu^C,. 

(5.70) 

Differentiating  Eqs.  (5.68-5.70)  produces 


dXc 

djs 


dXc 


dXc 


dX, 


__  =  C3(l)^  +  C,,(2)^  +  C,(:i)^ 

dx,  dx,  dx^.  dx, 

=  C.(l)^  +  C,(2)^  +  C,(3)-^'“ 


dXc  ’  dXc  dXc 

Lastly,  Eq.  (5.52)  can  be  rewritten 

Till  =  unit  ((r/rjvc  -  (r/t>,)rc) 
=  unit(C5) 


dXc 


which  can  be  differentiated  to  yield 

dTiLx  ^  ^ 

aXc  ViiCsIl  wcsrIdXc 

where 


dCi 

dX, 


^  T  *  w  *  *  'I'  ^  w  “*7' 

2Vcrc  -  r/ Vcl  -  r,v,'  r/r,/-r,.r,' 


Further,  differentiating  Eqs.  (.5..5;{)  and  (5.5-1)  yield 


dTw, 

(L 

fc  X  tic)(f^  X 

Vcf 

I 

—  "I 

dXc 

( 

ii(f,  X  t>c)l! 

3 

IK^c 

X 

Vc)I|/ 

r  0  i  . 

0  - 

■*C 

Vc  ' 

X 

1 

f  t»« 
ej 

'  o 

ic 

Zc 

0 

-Xc 

.  Vc  -if 

0 

-Vc 

Xc 

0  . 

dTiui 

dXc 

TcTJ  / 

llrdP  IlfcII 

; 

0 

(5.71) 

(5.7-2) 

(5.73) 

(5.74) 

(5.75) 

(5.76) 


(5.77) 


(5.78) 

(5.79) 


Tables  5.1  and  5.2  summarize  the  measurement  -  state  relationships. 
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Table  5.1;  Measurement  -  State  Relationships 


Parameter 

Derivative 

h\k 

Eq.  (5.40) 

dXc 

Eq. 

(5.42) 

R 

Eq.  (5.47) 

dR 

dX. 

Eq. 

(5.49) 

R 

Eq.  (5.48) 

dii 

dXc 

Eq. 

(5.50) 

Eq.  (5.60) 

dXc 

Eq. 

(5.62) 

d 

Eq.  (5.61) 

de 

dXc 

Eq. 

(5.63) 

5.2  Simulation  Results 

A  nonlinear  simulation  was  developed  to  model  the  rendezvous  nav¬ 
igation  filter  (see  Chapter  7).  Several  profiles  were  exeniled  in  order  to  assess 
the  performance  of  the  rendezvous  navigation  filter.  For  each  profile,  clock¬ 
wise  orbits  were  used.  A  planar  circular  orbit  was  also  used  to  insure  final 
position  errors  due  to  the  rendezvous  targeting  law  were  less  than  one  meter. 
The  incline  angle,  initial  condition  angle  and  transfer  time  were  fixed  at  the 
baseline  values  which  resulted  in  an  “in-plane”  rendezvous  profile. 

The  chaser  vehicle  was  initialized  19  km  behind  the  target  vehicle  in  a 
coelliptic  orbit.  An  error  vector  was  then  added  to  the  chaser  vehicle  to  model 
the  uncertainty  in  the  state  just  prior  to  the  77  maneuver.  Error  vectors  were 
generated  randomly  using  error  statistics  consistent  with  GIngiss  (35).  Thir¬ 
teen  particular  error  vectors  were  used  in  this  research:  each  one  corresponded 
to  a  constant  \a  error  from  the  6-dimensional  error  ellipsoid  at  the  77  luaneii 
ver.  Potter  [42]  and  Hitzl  [43]  both  describe  the  lechiiiqui*  u.s<‘(l  in  choosing 
these  vectors.  This  technique  consists  of  computing  a  scalar  using  the  error 
vector  and  the  inverse  of  the  covariance  matrix.  This  scalar  is  proportional 
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Table  5.2:  Parameters  used  by  Measurement  -  State  Relationship.s 


Parameter  Derivative 


R 

Eq.  (5.45) 

dR 

dXc 

Eq.  (5.56) 

V 

Eq.  (5.46) 

dV 

dX. 

Eq.  (5.57) 

a 

Eq.  (5.58) 

da 

dXc 

Eq.  (5.64) 

/? 

Eq.  (5.59) 

#- 

dXc 

Eq.  (5.65) 

*5 

Eq.  (5.68) 

du 

dX] 

Eq.  (5.71) 

h 

Eq.  (5.69) 

|k 

dX, 

Eq.  (5.72) 

its 

Eq.  (5.70) 

dk, 

dX. 

Eq.  (5.73) 

c. 

Eq.  (5.55) 

.N/A 

N/A 

Cj 

Eq.  (5.67) 

N/A 

N/A 

Ca 

Eq.  (5.67) 

N/A 

N/A 

c. 

Eq.  (5.67) 

N/A 

N/A 

c. 

Eq.  (5.75) 

c 

Eq.  (5.77) 

T/u 

Eq.  (5.52) 

aiv 

dX. 

Eq.  (5.76) 

TiL2 

Eq.  (5.53) 

dX. 

Eq.  (5.78) 

Til3 

Eq.  (5.54) 

dX. 

Eq.  (5.79) 

Figure  5.3:  Position  Error  Vectors  for  Monte  Carlo  .Analysis 

to  the  n-dimensional  probability  density  function  and  hence,  (onslanl  val¬ 
ues  of  this  scalar  equate  to  constant  surfaces  on  the  error  ellipsoid.  Finally, 
these  vectors  were  chosen  to  span  all  directions  of  the  77  error  ellipsoid.  Fig¬ 
ures  5.3  and  5.4  show  the  direction  of  these  position  and  velocity  error  vectors, 
respectively. 

The  proximity  sensor  wm  modeled  after  the  Space  Shuttle  rendezvous 
radar  [44].  Sensor  measurements  consisted  of  relative  range  and  line-of-sight 
angles  and  were  taken  every  minute.  Measurement  statistics  and  sensor  capa¬ 
bilities  were  chosen  similar  to  the  Space  Shuttle  rendezvous  radar  [44]  and  were 
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consistent  with  an  advanced  rendezvous  tracking  system  study  [45].  Random 
noise  was  added  to  each  measurement,  but  measurement  biases  and  unmod¬ 
eled  accelerations  were  omitted  from  the  estimation  v  .,r. 

Figure  5.5  shows  the  average  state  estimation  error  of  all  13  error 
vectors  (solid  line)  in  the  LVIR  coordinate  frame;  also  shown  is  the  standard 
deviation  of  the  state  (dotted  lines).  Appendix  C.l  contains  ihe.sc  same  data 
for  each  individual  simulation.  State  estimation  errors  converged  quickly  after 
only  a  few  measurements  were  incorporated. 

In  Chapter  4,  midcourse  correction  maneuvers  were  shown  to  sub¬ 
stantially  reduce  final  position  errors  due  to  the  rendezvous  targeting  law;  they 
also  aid  in  reducing  the  effects  of  navigation  errors.  The  optimal  placement  of 
a  midcourse  correction  maneuver  within  the  rendezvous  profile  depends  on  the 
navigation  errors  at  the  time  of  the  maneuver  and  the  sensitivity  of  the  final 
position  errors  at  the  time  of  rendezvous  with  respect  to  maneuver  velocity 
errors.  Two  parametric  studies  were  completed  to  characterize  the  optimal 
placement  of  a  midcourse  correction  maneuver. 

5.2.1  Navigation  Errors  Parametric  Study 

The  first  study  investigated  was  the  navigation  errors  along  the  ref¬ 
erence  trajectory  for  the  baseline  “in-plaiie”  rendezvous  profile.  This  study 
was  designed  to  answer  the  question,  “Is  there  a  place  in  the  rendezvous  pro¬ 
file  where  the  navigation  errors  are  small  in  the  direction  of  the  AV'  required 
from  a  midcourse  correction  maneuver?”  If  such  a  place  exists,  a  midcourse 
correction  maneuver  should  be  executed  at  this  point  in  the  rendezvous  profile 
because  the  midcourse  correction  maneuver  is  insensitive  to  navigation  errors. 

To  determine  if  such  a  position  in  the  rendezvous  profile  exists,  a 


covariance  analysis  was  done.  Figure  5.6  (a)  shows  the  time  history  of  the 
magnitude  of  the  principle  axes  of  the  LVIR  covariance  position  error  ellipsoid. 
Also  indiciated  on  Fig.  5.6  (a)  is  the  general  correlation  between  the  LVIR 
coordinate  frame  and  the  principle  axes.  Figure  5.6  (b)  gives  the  LVIR  azimuth 
and  elevation  angles  of  the  largest  principle  axis  of  the  covariance  position 
error  ellipsoid.  Figure  5.7  provides  these  same  data  for  the  LVIR  covariance 
velocity  error  ellipsoid. 

The  covariance  position  error  ellipsoid  is  initialized  as  a  sphere.  With 
the  first  couple  of  proximity  sensor  measurements,  the  sphere  shrinks  at  an 
essentially  uniform  rate.  Then,  for  a  brief  time  the  out-of-plane  direction 
is  the  largest  principle  axis,  but  then  it  decreases  and  the  radial  direction 
becomes  the  largest  principle  axis.  It  then  shrinks  and,  finally,  the  downtrack 
direction  becomes  the  largest  principle  axis.  The  final  covariance  position 
error  ellipsoid  is  aligned  with  the  LVIR  coordinate  frame,  but  rotated  4.9  deg 
about  the  positive  out-of-plane  axis. 

The  covariance  velocity  error  ellipsoid  behaves  in  an  similar  manner 
except  for  the  out-of-planc  direction  never  being  the  largest  principle  axis.  The 
final  covariance  velocity  error  ellipsoid  is  also  aligned  with  the  LVIR  coordinate 
frame,  but  it  is  rotated  5.0  deg  about  the  positive  out-of-plane  axis. 

Note  that  even  though  these  error  ellipsoids  lose  their  spherical  shape, 
which  might  suggest  a  desirable  location  for  a  iiiidcourse  cuireclioii  iiiaiieiiv«i , 
the  absolute  magnitude  of  all  of  the  principle  axes  is  small  after  only  a  few 
measurements  have  been  incorporated  into  the  rendezvous  navigation  filter. 
Thus,  placement  of  a  midcourse  correction  maneuver  is  not  directly  dependent 
on  navigation  errors. 


Angle  (deg)  ^  Velocity  (m/s) 


13 


(b)  Angie  to 


Axis 
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5.2.2  Final  Position  Error  Sensitivity  Parametric  Study 


The  second  study  investigated  was  the  sensitivity  of  the  final  position 
errors  with  respect  to  maneuver  velocity  errors.  This  study  was  designed  to 
answer  the  question,  “Is  there  a  place  in  the  rendezvous  profile  where  final 
position  errors  are  insensitive  to  velocity  errors  from  executing  a  midcourse 
correction  maneuver?”  If  such  a  place  exists,  a  midcourse  correction  maneuver 
should  be  executed  at  this  point  in  the  rendezvous  profile  because  any  velocity 
errors  in  the  maneuver  execution  will  have  less  affect  on  final  position  errors. 


To  determine  if  such  a  position  in  the  rendezvous  profile  exists,  a 
sensitivity  matrix  was  computed;  the  sensitivity  matrix  is  defined  as 


dSpTF  A 
dV 


r  ^  ^  dSj 
W,  W,  W, 

^  ^ 

Wi  Wr 

d6r  dSr  ^ 


(5.80) 


where  the  (  )j,  (  )„  and  (  ),  subscripts  indicate  the  dowiitrack.  oul-of- 

plane  and  radial  directions,  respectively.  Each  partial  derivative  in  Eq.  (5.80) 
was  computed  numerically  along  the  reference  trajectory  for  the  baseline  “in¬ 
plane”  rendezvous  profile  using  the  central  differences  technique  with  four 
evaluations. 


Figure  5.8  gives  a  time  history  of  the  largest  singular  value  of  the 
sensitivity  matrix;  individual  time  histories  of  the  elements  of  the  sensitivity 
matrix  are  included  in  Appendix  C.2.  Figure  5.8  shows  no  significant  sen¬ 
sitivity  of  the  final  position  errors  with  respect  to  velocity  errors.  However, 
it  does  substantiate  the  general  philosophy  that  the  final  position  errors  will 
be  smaller  given  a  smaller  transfer  time  between  the  niidcour.se  correction 
maneuver  and  the  time  of  rendezvous.  Nevertheless,  th»*  placement  of  a  mid 
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course  correction  maneuver  is  not  directly  dependent  on  final  position  error 
sensitivities. 

Thus,  since  the  placement  of  a  midcourse  correction  maneuver  is  not 
directly  dependent  on  navigation  errors  or  final  position  error  sensitivities,  the 
midcourse  correction  maneuvers  can  be  placed  in  the  rendezvous  profile  based 


Chapter  6 


Applications 


To  demonstrate  the  entire  guidance  and  navigation  system,  this 
chapter  considers  two  profiles  for  terminal  phcise  rendezvous  in  the  planar 
circular  orbit  and  two  equivalent  profiles  in  a  'H-i  guided  halo  orbit.  The  two 
profiles  chosen  correspond  to  “in-plane”  (profile  1)  and  “out-of-plane”  (profile 
2)  relative  trajectories  for  the  planar  circular  orbit;  the  initial  condition  angles 
were  46.08  deg  and  180.00  deg,  respectively.  These  same  initial  condition  an¬ 
gles  were  used  in  the  H-i  guided  halo  orbit  demonstration.  The  transfer  time 
was  also  held  constant  at  9  hrs  for  both  demonstrations.  The  incline  angle  for 
the  planar  circular  orbit  demonstration  was  22.85  deg;  this  is  the  same  as  the 
incline  angle  of  the  'Hj  guided  halo  orbit  used  in  this  research. 

Monte  carlo  simulations  >  ere  done  using  the  13  error  vectors  se¬ 
lected  in  Chapter  5.  For  comparison  purposes,  a  reference  simulation  was 
also  done  assuming  “ideal”  navigation.  Four  midcourse  correction  maneuvers 
were  added  to  the  rendezvous  profile  to  correct  for  rendezvous  targeting  and 
rendezvous  navigation  filter  state  estimation  errors.  These  maneuvers  were 
placed  on  time  at  TI  plus  two,  Tl  plus  four,  Tl  plus  six  and  Tl  plus  eight 
hours.  The  Space  Shuttle  program  also  uses  four  midcourse  correction  ma¬ 
neuvers  in  its  rendezvous  profile;  three  of  these  maneuvers  are  placed  on  lime 
while  one  is  executed  on  an  elevation  angle  to  the  target  vehicle. 

This  chapter  presents  a  summary  of  these  simulations;  Appendix  D 
tabulates  the  propulsion  requirements  and  final  miss  distance  for  each  indi- 
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vidual  rendezvous  profile.  Propulsion  requirements  are  given  in  terms  of  the 
velocity  required  for  each  maneuver.  Final  miss  distance  is  defined  as  the  error 
from  the  desired  position  offset  at  the  time  of  rendezvous. 


6.1  Planar  Circular  Orbit  Demonstration 


Figure  6.1  shows  the  reference  relative  trajectory  for  the  “in-plane” 
profile;  Fig.  6.2  applies  to  the  “out-of-plane”  profile.  In  both  profiles,  the 
chaser  vehicle  initiated  the  rendezvous  profile  19  km  behind  the  target  vehicle 
in  a  coelliptic  orbit.  These  relative  trajectories  have  two  main  differences. 
In  profile  1,  the  “in-plane”  relative  trajectory  has  a  maximum  out-of-plane 
position  less  than  20  m  and  goes  only  slightly  below  the  target  vc'liicle.  This 
relative  trajectory  is  essentially  straight  line  motion  and  shows  little  orbital 
mechanics  effects.  By  contrast,  the  maximum  out-of-plane  relative  position  in 
profile  2  is  greater  than  800  m  and  the  “out-of-plane”  relative  trajectory  goes 
below  the  target  vehicle  over  6  times  the  amount  of  the  “in-plane”  relative 
trajectory.  Orbital  mechanics  effects  are  clearly  present  in  this  rendezvous 
profile. 


These  rMative  trajectory  differences  are  also  reflected  in  the  propul¬ 
sion  requirements  and  final  miss  distance  of  each  rendezvous  profile.  Table  6.1 
summarizes  the  propulsion  requirements  for  both  the  planar  circular  orbit  and 
the  7^2  guided  halo  orbit;  the  reaction  control  system  (RCS)  fuel  consumed 
was  computed  using 


Wp  = 


(6.1) 


The  initial  spacecraft  weight  was  assumed  to  be  the  maximum  allowed  (13,000 
lbs)  under  the  assumptions  of  the  restricted  three-body  problem.  Also,  a 
medium  performance  RCS  jet  (/,p  =  270  sec)  was  assumed;  this  is  a  typi- 


Figure  6.1:  Planar  Circular  Orbit  Reference  Relative  Trajectory 

(IC=46.08  deg,  INC=22.85  deg,  TT=9  hrs) 
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Table  6.1:  Propulsion  Requirements  Summary 


Error 

Rendezvous  Maneuver  (m/s) 

RCS 

Vector 

TI 

MCCl 

MCC2 

MCCi 

MCCi 

TF 

Total 

Fuel  (lbs) 

Circle  #1 
Reference 
Average 

0.686 

0.686 

0.000 

0.041 

0.000 

0.010 

0.000 

0.010 

0.000 

0.014 

0.491 

0.493 

1.177 

1.254 

5.776 

6.153 

Circle  #2 
Reference 

0.680 

0.000 

0.000 

0.000 

0.000 

0.643 

1.322 

6.487 

Average 

0.680 

0.041 

0.009 

0.008 

0.014 

0.634 

1.384 

6.791 

Halo#l 

Reference 

0.627 

0.019 

0.022 

0.027 

0.044 

0.684 

1.421 

6.972 

Average 

0.627 

0.041 

0.023 

0.032 

0.042 

0.677 

1.441 

7.071 

Halo  #2 
Reference 

0.704 

0.008 

0.008 

0.009 

0.013 

0.672 

1.413 

6.933 

Average 

0.704 

0.042 

0.013 

0.012 

0.021 

0.663 

1.455 

7.139 

cal  Space  Shuttle  RCS  jet.  Table  6.2  shows  the  final  miss  distance  for  both 
demonstrations.  Velocity  requirements  increased  12.3%  from  the  “in-plane” 
relative  trajectory  to  the  “out-of- plane”  relative  trajectory:  final  nnss  distance 
increased  21.3%. 

State  estimation  errors  from  the  rendezvous  navigation  filter  were- 
also  small.  Figures  6.3  and  6.4  show  the  average  state  estimation  error  of  all 
13  error  vectors  (solid  line)  and  the  associated  1<7  standard  deviation  of  the 
state  (dotted  line).  These  data  are  given  in  the  LVIR  coordinate  frame.  The 
rendezvous  navigation  filter  was  not  substantially  affected  by  the  geometry 
of  the  relative  trajectory.  These  state  estimation  errors  caused  the  velocity 
requirements  to  increase  6.5%  and  4.7%  for  the  “in-plane”  and  “out-of-plane” 
profiles,  respectively.  Final  miss  distance  also  increased  to  12.7  m  and  15.3  m 
for  these  profiles,  respectively. 

Even  though  these  increases  due  to  relative  t*  ajectory  differences  and 
rendezvous  navigation  filter  state  estiuiatiuii  errors  are  significant,  it  should 
be  noted  that  the  absolute  propulsion  re<|uireinents  anil  final  miss  clistame 


Figure  6.3:  Planar  Circular  Orbit  Average  Chaser  Vehicle  State  Estimation 
Error  (IC=46.08  deg,  INC=22.85  deg,  TT=9  hrs) 


Figure  6.4:  Pl&nar  Circular  Orbit  Average  Chaser  Vehicle  State  Estimation 
Error  (IC=180.00  deg,  1NC=22.85  deg,  TT=9  hrs) 
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Table  6.2:  Final  Miss  Distance  Summary 


Final  Miss  Distance  (m) 

Downtrack  Out-of- Plane 

Radial 

Total 

Circle  #1 

Reference 
Average  Magnitude 

0.000 

4.919 

-0.001 

7.204 

0.000 

6.983 

0.001 

12.666 

Circle  #2 

Reference 

0.000 

0.000 

-0.001 

0.001 

Average  Magnitude 

8.086 

9.779 

5.686 

15.362 

Halo  #1 

Reference 

-7.67.3 

19.394 

0.124 

20.857 

Average  Magnitude 

9.160 

19.585 

9..303 

2.5.924 

Halo  #2 

Reference 

-4.834 

-3.191 

-0.885 

5.859 

Average  Magnitude 

5.898 

10.222 

5.036 

14.109 

were  small  in  each  profile. 


6.2  7i2  Guided  Halo  Orbit  Demonstration 

Figure  6.5  shows  the  reference  relative  trajectory  for  profile  1; 
Fig.  6.6  applies  to  profile  2.  In  both  profiles,  the  chaser  vehicle  initiated 
the  rendezvous  profile  approximately  19  km  behind  and  slightly  above  the 
target  vehicle.  The  trajectories  clearly  show  the  same  characteristic  as  the 
planar  circular  orbit  with  respect  to  the  radial  component  of  the  rendezvous 
profile.  Not-so-obvious  is  why  the  maxiniuin  out  uf-plane  position  diil  not 
follow  the  same  trend.  This  exception  is  explained  by  noting  that  for  the 
planar  circular  orbit,  the  (-axis  velocity  component  of  the  chaser  vehicle  was 
zero  prior  to  the  T I  maneuver  and  after  the  77  maneuver,  the  chaser  vehicle 
moved  in  the  negative  (  direction.  However,  for  the  Tij  guided  halo  orbit, 
the  {-axis  velocity  component  of  the  chaser  vehicle  was  positive  prior  to  the 
T I  maneuver.  After  the  T I  maneuver,  the  chaser  vehicle  again  moved  in  the 


(d)  ReiativeTtajeciory 


Figure  6.5:  Wj  Guided  Halo  Orbit  Reference  Relative  Trajectory 
(IC=46.08  deg,  TT=9  hrs) 
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negative  (  direction.  Hence,  the  pre-T I  ^-axis  velocity  opposed  a  portion  of 
the  post-T/  motion  and  less  total  out-of-plane  relative  motion  resulted.  In 
both  profiles,  orbital  mechanics  effects  are  seen. 

The  effect  of  these  relative  trajectory  differences  was  also  present 
in  the  propulsion  requirements  and  final  miss  distance.  Essentially  the  same 
trends  as  for  the  planar  circular  orbit  demonstration  were  seen.  State  estima¬ 
tion  errors  from  the  rendezvous  navigation  filter  were  also  small  just  as  in  the 
planar  circular  orbit  demonstration;  Figs.  6.7  and  6.8  give  the  average  state 
estimation  error  of  all  13  error  vectors  (solid  line)  and  the  associated  l<r  stan¬ 
dard  deviation  of  the  state  (dotted  line).  These  data  are  given  in  the  LVIR 
coordinate  frame.  Again,  the  rendezvous  navigation  filter  performance  was 
not  affected  by  the  geomc'ry  of  the  relative  trajectory.  The  state  estimation 
errors  caused  the  velocity  requirements  to  increase  1.4%  and  3.0%  for  profiles 
1  and  2,  respectively.  Final  miss  distance  also  increased  24.3%  and  140.8% 
for  these  profiles,  respectively.  When  considering  the  absolute  magnitude  of 
the  parameters,  these  increases  are  approximately  the  saiiu'  as  for  the  planar 
circular  orbit  demonstration. 

As  before,  these  increases  are  sigiiifica^'t,  but  it  should  be  emphasized 
that  the  absolute  propulsion  requirements  and  final  miss  distance  remained 
small  in  both  profiles. 


20. 


Figure  6.7:  'H2  Guided  Halo  Orbit  .\verage  Chaser  Vehicle  Stale  Estimation 
Error  {IC=46.08  deg,  TT=9  hrs) 


Chapter  7 


Nonlinear  Simulation  Package 


This  simulation  package  consists  of  three  executable  executive  pro¬ 
grams,  associated  subroutines,  and  numerous  data  plotting  command  files. 
The  pre-processor  program  is  an  interactive  parameter  input  program  used  to 
define  approximately  500  selectable  parameters  required  by  the  main  simula¬ 
tion.  The  post-processor  program  edits  stored  data  and  computes  variables 
of  interest  such  as  the  statistics  on  the  measurement  residuals  from  the  ren¬ 
dezvous  navigation  filter.  All  source  code  is  written  in  FORTRAN;  the  data 
plotting  command  files  use  the  XGRAMI  application. 

This  chapter  begins  by  describing  the  general  capabilities  of  the  sim¬ 
ulation  package.  The  functional  modules  of  the  main  simulation  are  then 
discussed.  Lastly,  the  pre-processor  and  post- processor  programs  are  outlined. 

7.1  General  Capabilities 

The  scenario  being  simulated  is  the  terminal  phase  rendezvous  be¬ 
tween  a  chaser  and  target  vehicle  in  the  circular  restricted  three- body  problem. 
The  target  and  chaser  vehicles  are  assumed,  initially,  to  be  in  small  radius  halo 
orbits  about  the  translunar  equilibrium  point  in  the  Earth-Moon  three-body 
system.  Parameters  identifying  the  specific  three-body  system  are  selectable. 
The  type  of  translunar  halo  orbit  is  selectable  between  a  planar  circular  orbit 
and  a  W2  guided  halo  orbit.  If  a  Kj  guided  halo  orbit  is  selected,  the  target 
and  chaser  vehicles  can  be  placed  in  independent  orbits.  If  a  planar  circular 
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orbit  is  selected,  the  chaser  vehicle  is  initialized  at  a  selectable  offset  position 
from  the  target  vehicle  in  a  coelliptic  orbit. 

The  main  simulation  computes  TI  and  TF  maneuvers  plus  up  to 
eight  midcourse  correction  {MCCl  —  MCC8)  maneuvers.  All  maneuvers  are 
modeled  impulsively.  The  TI  maneuver  may  occur  any  lime  after  llie  inilial 
state  vector  time  for  the  chaser  vehicle.  Up  to  five  distinct  sensor  passes  may 
occur  prior  to  each  maneuver  except  T /.  Times  for  all  maneuvers  plus  sensor 
pass  start  and  stop  times  are  selectable.  The  initial  vector  state  vector  time 
for  the  chaser  and  target  vehicles  and  the  TI  time-of-ignition  (TIG)  are  in 
absolute  time;  all  other  times  are  relative  to  TI  TIG. 

Five  coordinate  frames  are  used  in  the  simulation;  all  computations 
are  done  in  the  inertial  coordinate  frame. 

1.  Inertial  coordinate  frame:  Centered  at  the  translunar  equilibrium  point 
with  the  ^-axis  pointing  along  a  line  from  the  Earth  to  the  Moon,  the  tj- 
axis  in  the  Earth- Moon  plane  and  pointing  in  the  direction  ol  tlu'  Moon  s 
orbital  motion,  and  the  (.'-axis  completing  the  right  handed  (  ounliiiale 
system. 

2.  Local  vertical  inertial  rectangular  (LVIR)  coordinate  frame;  Centered 
at  either  vehicle  with  the  radial-axis  along  the  negative  position  vec¬ 
tor,  the  out-of-plane-axis  along  the  crossproduct  vector  of  the  negative 
position  and  velocity  vectors,  and  the  downtrack-axis  completing  the 
right-handed  coordinate  system. 

3.  Target  centered  curvilinear  coordinate  frame:  Centered  at  the  target  ve¬ 
hicle.  The  radial  distance  is  defined  as  the  distance  between  the  target 
vehicle,  when  propagated  in  a  circular  orbit  to  a  matching  |)ha.s<“  aiigh- 
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with  the  chaser  vehicle,  and  the  chaser  vehicle  in  a  radial  direction.  The 
out-of-plane  distance  is  defined  as  the  distance  between  the  target  vehi¬ 
cle  and  the  chaser  vehicle  in  the  out-of-plane  direction.  The  downtrack 
distance  is  defined  as  the  curvilinear  distance  along  the  propagated  cir¬ 
cular  target  orbit.  The  positive  direction  for  all  three  axes  is  defined  in 
the  same  manner  as  the  LVIR  coordinate  frame. 

4.  Body  centered  coordinate  frame:  Centered  at  the  chaser  vehicle  and 
transformed  from  the  LVIR  coordinate  frame  by  a  yaw-pitch-roll  set  of 
Euler  angles. 

5.  Sensor  centered  coordinate  frame:  Centered  at  a  specified  point  in  the 
body  centered  coordinated  frame  and  transformed  from  the  body  cen¬ 
tered  coordinate  frame  by  a  yaw-pitch-roll  set  of  Euler  Angles. 

7.2  Main  Simulation  Functional  Description 

The  main  simulation  consists  of  five  primary  modules  plus  several 
support  modules.  Figure  7.1  shows  an  overview  of  how  these  five  modules 
interact  to  form  the  guidance  and  navigation  system.  Figures  7.2  and  7.3 
summarize  how  the  target,  chaser  and  estimated  trajectories  are  initialized  and 
integrate  from  the  initial  condition  {IC)  to  the  beginning  of  tlx-  rendc'zvons 
profile  (TJ)  and  finally  to  the  time  of  rendezvous  (77* ). 

7.2.1  Target  Vehicle  Environment  Module 

Figure  7.4  gives  the  functional  description  of  the  target  vehicle  mod¬ 
ule.  The  target  vehicle  can  be  placed  in  an  guided  halo  orbit  or  in  a  planar 
circular  orbit  by  a  selectable  flag.  In  either  case,  the  target  vehicle  remains 
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Figure  7.2:  Guided  Halo  Orbit  Trajectories 
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Figure  7.3:  Planar  Circular  Orbit  Trajectories 
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(a)  Target  Vehicle  in  Guided  Halo  Orbit 


Figure  7.4:  Target  Vehicle  Environment  Functional  Description 
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passive  with  respect  to  the  rendezvous  profile. 

For  a  'H2  guided  halo  orbit,  the  target  vehicle’s  initial  states  and 
time  are  selectable.  Orbital  rotation  is  dependent  on  these  states.  The  target 
vehicle  dynamics  consists  of  the  circular  restricted  three-body  equations  of 
motion.  The  amount  of  disturbances  added  to  the  target  vehicle  dynamics  is 
selectable.  Position  measurements  for  the  guidance  law  are  generated  assum¬ 
ing  a  “ideal”  guidance  law  sensor.  The  amount  of  measurement  noise  added 
to  these  “ideal”  measurements  is  selectable.  The  fixed  power  spectrum  for  the 
disturbances  and  measurement  noise  is  selectable  but  must  be  consistent  with 
the  state  space  realization  of  the  guidance  law.  Further,  the  non-dimensional 
frequency  and  radius  of  the  deseed  halo  orbit  are  selectable  but  must  also  be 
consistent  with  the  guidance  law. 

If  a  planar  circular  orbit  is  selected,  the  orbit’s  orientation  with  re¬ 
spect  to  the  r}-(l  plane  (incline  angle)  is  specified  by  a  selectable  distance.  Its 
orbital  rotation  (clockwise  or  counterclockwise)  is  controlled  through  a  sepa¬ 
rate  selectable  parameter.  The  initial  states  are  computed  automatically  using 
the  selectable  non-dimensional  frequency  and  radius  of  tin'  halo  orbit  in  addi 
tion  to  its  initial  time.  The  initial  time  of  the  target  vehicle  is  (  omputed  using 
the  selectable  cha.'ver  vehicle’s  77  downtrack  offset  position  and  the  target 
vehicle’s  non-dimensional  frequency  and  halo  orbit  radius,  rin*  target  vehi¬ 
cle’s  dynamics  consists  of  simple  second  order  linear  time  invariant  differential 
equations  that  are  only  a  function  of  the  non-dimensional  freijuency. 

7.2.2  Chaser  Vehicle  Environment  Module 

Figure  7.5  depicts  the  functional  description  of  the  chaser  vehicle 
module  during  the  rendezvous  profile;  Fig.  7.4  (a)  would  apply  to  the  chaser 
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(b)  Target  Vehicle  in  Planar  Circular  Orbit 

Figure  7.5:  (’ha.scr  Vehicle  Kuviroiiuieiit  Fuiiclioiial  i)(‘S(  ri|>lu)ii 

vehicle  prior  to  TI  if  the  target  vehicle  is  in  a  guided  halo  orbit.  In  this 
case,  the  chaser  vehicle  is  initialize  in  an  independent  guided  halo  orbit 
(sec  Fig.  7.2).  The  initial  chaser  vehicle  state  and  time  are  selectable.  For 
integration  from  the  initial  time  to  77,  the  chaser  vehicle  dynamics  consists  of 
the  circular  restricted  three-body  equations  of  motion  plus  a  Hj  optimal  guid¬ 
ance  law.  The  amount  of  disturbances  added  to  the  cha-ser  vehicle  dynamics 
is  selectable.  The  position  measurements  for  the  guidance  law  are  generated 
assuming  a  “ideal”  guidance  law  scn.sor.  The  amount  of  measurement  iioi.se 
added  to  these  “ideal”  measurements  is  .s<*h*ctable.  The  li.xe<l  power  spi*< 
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trum  for  the  disturbances  and  measurement  noise  is  selectable  but  must  be 
consistent  with  the  state  space  realization  of  the  guidance  law.  Further,  the 
non-dimensional  frequency  and  radius  of  the  desired  halo  orbit  are  selectable 
but  must  also  be  consistent  with  the  guidance  law. 

During  the  rendezvous  profile,  the  chaser  vehicle  is  assumed  to  be 
in  coasting  flight.  The  initial  state  and  time  are  taken  as  the  last  state  prior 
toTI  plus  a  table  look-up  error  vector  and  TI  TIG,  respectively.  The  chziser 
vehicle  dynamics  consists  of  the  circular  three-body  equations  of  motion.  The 
amount  of  disturbances  added  to  the  chaser  vehicle  dynamics  is  selectable 
and  is  independent  of  the  pre-TI  phase.  The  fixed  power  spectrum  for  the 
disturbances  is  also  selectable  and  independent  of  the  prc-  I'l  phase. 

If  the  target  vehicle  is  in  a  planar  circular  halo  orbit,  the  <  baser 
vehicle  is  initialize  at  T/  based  on  selectable  offsets  from  the  target  vehicle 
assuming  a  coelliptic  orbit.  A  table  look-up  error  vector  is  added  to  this 
computation  (see  Fig.  7.3).  The  chaser  vehicle’s  time  is  assumed  to  be  TI  TIG. 
The  chaser  vehicle  is  assumed  to  be  in  coasting  flight,  its  dynamics  consists 
of  either  the  nonlinear  or  linear  circular  restricted  three-body  equations  of 
motion  depending  on  a  selectable  flag. 

7.2.3  Proximity  Sensor  Module 

Figure  7.6  shows  tin*  functional  rlescriptimi  of  ihe  prtrximily  sensor 
module.  The  proximity  sensor  module  is  niodel(>d  after  I  In-  Space  Shuttle  ten 
dezvous  radar  [l  l]  and  considers  four  types  of  ineasnrenient.s;  relative  laiige. 
relative  range  rate,  and  two  line-uf-sighi  angh*s.  Tin'  hn«‘-ol’  sight  angles  used 
by  the  rendezvous  navigation  filter  an*  roll  and  pitch  angles  in  the  body  cen 
tercd  coordinate  frame  and  are  computed  <lin“<  tly  from  tin  shaft  and  trunnion 


Figure  7.6;  Proximity  Sensor  Functional  Description 
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angles  in  the  sensor  centered  coordinate  frame. 

All  measurements  are  generated  regardless  of  whether  or  not  they 
are  used  by  the  rendezvous  navigation  filter.  The  “ideal”  measurements  are 
checked  by  the  sensor  capabilities  model  to  ensure  that  they  are  within  sensor 
hardware  limits;  these  limits  are  selectable.  If  the  generated  measurement  is 
outside  of  the  sensor  hardware  limits,  the  data  good  flag  a.ssociated  with  that 
measurement  type  is  set  to  “bad.” 

The  addition  of  random  noise  to  the  “ideal”  ineasureinents  is  duiu' 
in  the  sensor  error  model  and  is  selectable.  Characteristics  of  the  random 
noise  are  also  selectable.  The  sensor  error  model  also  has  the  capability  to 
add  random  biases  to  each  measurement.  This  addition  is  independent  of 
the  random  noise  addition  and  is  selectable.  The  random  bias  is  computed 
only  once  per  phase;  phases  are  defined  by  the  range  to  the  target  vehicle. 
Characteristics  of  the  random  bia.ses  are  also  selectable. 

7.2.4  Rendezvous  Navigation  Module 

Figure  7.7  gives  the  functional  description  of  the  r<'n(lezvous  naviga¬ 
tion  filter.  The  rendezvous  navigation  module  is  formulated  as  a  < out  imioiis 
dynamics,  discrete  measurements  e.xt<*nded  Kalman  lilt(“i  .  .Measurements  are 
generated  by  the  proximity  sensor  module  ami  may  be  im  orporated  into  the 
filter  in  any  combination;  that  is,  ndative  range  and  line-of-sight  angles,  or 
relative  range  and  range  rate,  etc.  The  rendezvous  navigation  filter  clnrks  the 
data  good  flag  of  each  measurement  used;  if  any  data  good  flag  of  a  measure¬ 
ment  selected  is  “bad,”  the  entire  s<*t  of  incasurem<‘nts  is  rejerted. 

The  estimated  states  consist  of  the  clia,ser  vehicle  iin'rtial  state  plus 
either,  both  or  neither  of  unmodeled  accelerations  acting  on  the  chaser  vehi- 
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cle  and  measurement  biases.  Which  states  are  estimated  is  selectable.  I'he 
unmodeled  accelerations  and  measurement  biases  are  each  modeled  as  a  first 
order  Gauss-Markov  process  with  selectable  time  constants.  Initialization  of 
the  estimated  chaser  vehicle  states  is  either  the  pre-Tl  chaser  vehicle  states 
for  a  H2  guided  halo  orbit  or  the  pre-TI  target  vehicle  states  plus  an  offset 
position  for  a  planar  circular  orbit  (see  Figs.  7.2  and  7.3).  Initial  values  of  the 
unmodeled  accelerations  and  measurement  biases  are  selectable.  The  initial 
time  associated  with  the  estimation  states  is  T/  TIG. 

The  initial  covariance  ntatrix  assumes  a  diagonal  form  and  has  a 
selectable  standard  deviation  for  each  estimated  state.  The  values  associated 
with  the  estimated  chaser  vehicle  states  are  used  to  generate  the  table  look¬ 
up  error  vectors  for  initialization  of  the  chaser  vehicle  environiiKMit  I'he  fixed 
power  spectrum  matrix  for  the  proct^s  and  measur(*m<'iit  noise  ti'iins  assume 
a  diagonal  form  and  have  a  selectable  standard  <h>viation  for  <‘a(  h  parametei. 

Each  sensor  pass  is  specified  by  start  and  stoj)  times,  nmiibei  of 
measurements  taken,  and  measurement  set  type.  All  of  these  parameters 
are  selectable  for  each  sensor  pass.  For  each  midcourse  corr«*ction  maneuver, 
the  covariance  matrix  is  enlarged  to  reflect  the  uncertainty  in  the  maneuver 
execution.  During  this  enlargement,  the  uncertainty  in  the  estimated  chaser 
vehicle  velocity  states  is  increased  by  a  selectalde  amount  while  the  overall 
correlation  coefficient  matrix  is  preserved. 

7.2.5  Rendezvous  Targeting  Module 

Figure  7.8  gives  the  functionaf  <h*scri[)tion  of  the  rendezvous  target 
ing  module.  The  7’/  ami  M('('  maneuvers  are  (oniputed  from  the  estimated 
chaser  vehicle  state  and  target  vehicle  stale  to  satisfy  the  <lesire<l  position  off 
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set  at  the  time  of  rendezvous.  The  position  offset  at  th('  tinu'  ol  reiide/.xous 
is  selectable. 

Two  rendezvous  targeting  laws  have  been  implemented  in  the  simu¬ 
lation  for  the  TI  and  MCC  maneuvers.  The  first  law  is  developed  in  Chapter 
4;  the  second  law  is  an  approximation  of  the  first  and  consists  of  straight 
line  motion.  Choosing  which  targeting  law  is  used  is  accomplished  with  a  se¬ 
lectable  flag.  An  iteration  cycle  eliminates  the  error  due  to  the  chaser  vehicle 
linearization  process  in  the  targeting  law  development;  the  iteration  tolerance 
and  maximum  number  of  iteration  cycles  are  selectable. 

The  TF  maneuver  is  computed  to  satisfy  the  desired  velocity  offsets 
at  the  time  of  rendezvotis;  these  offsets  are  selectable,  riu'  chHS('r  v('hi(  le  and 
target  vehicle  states  are  used  in  this  computation. 

7.2.6  Support  Modules 

A  Runge-Kutta.  7(8)  numerical  integrator  is  use<l  for  all  integration. 

The  initial  time  step,  absolute  error  and  relative  error  are  selectable.  Each 
integration  interval  may  be  broken  up  into  sub-intervals  for  data  storage.  The 
number  of  sub- intervals  is  selectable  and  applies  to  all  integration  intervals 
except  those  intervals  in  the  middle  of  a  sensor  pass. 

All  noises  used  in  this  simulation  are  g<*nerated  by  an  IMSL  pseudo¬ 
random  number  generator  from  a  staiulard  normal  distribution  using  an  acceptance- 
rejection  method.  The  seed  for  t!u'  l.MSb  psemio-raiidom  mmii)er  geiieratoi 
is  selectable. 

Available  output  from  the  main  simulation  toiisists  of  digital  data 
displays,  plots,  and  data  fih's.  Tin*  first  two<ligital  <latH  displays  provide  a  ma 
mmver  summary  of  the  rendezvous  profile;  each  maneuver  is  detaih‘<i  in  terms 
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of  the  velocity  required  and  relative  position  at  the  time  of  the  maneuver. 
Two  additional  displays  are  available  which  list  any  proximity  sensor  mea¬ 
surements  that  were  rejected  and  the  final  measurement  residual  statistics. 
Seventy  plots  are  available  which  show  the  relative  trajectory;  environment 
and  estimated  state  time  histories;  and  estimation  error,  covariance,  and  cor¬ 
relation  coefficient  time  histories.  Eighteen  raw  data  files  can  be  store  by  the 
main  simulation;  selectable  data  storage  flags  are  used  to  designate  which  data 
is  stored  for  each  individual  profile  executed. 

7.3  Pre-Processor  Program 

The  pre-processor  program  is  an  interactive  parameter  input  pro¬ 
gram  used  to  define  approximately  500  selectable  paramcMc'rs  rc'cjiiired  by  the 
main  simulation.  For  example,  these  parameters  sj)ecify  how  many  ic'iidcvvoiis 
maneuvers  and  sensor  passes  are  to  be  executed  as  well  as  the*  t  imc's  assoc  ialc‘cl 
with  these  maneuvers  and  passe's.  I'hey  ele*fine*  the*  initial  chaser  and  target 
vehicle  conditions  and  orbits.  They  also  control  the  capabilitie's  and  c  liaracter- 
istics  of  the  rendezvous  navigation  filter  and  the  proximity  se'nsor.  Table  7.1 
lists  the  main  menu  options  for  the  pre-proce'ssor  program.  When  the  pre¬ 
processor  program  is  exited,  two  data  files  are  stored  with  the  upelated  values 
of  all  parameters;  these  data  files  are  used  to  initialize  the  main  simulation. 

7.4  Post-Processor  Program 

The  post- processor  program  edits  ston'd  data  and  cc)iii|)ntes  vari¬ 
ables  of  interest  such  as  the'  statistics  on  tin*  iiH'asmc'ineni  rc'sidnals  from  the 
rendezvous  navigation  filter;  this  program  is  also  c  cnnplelelv  inlerac  l  i\c'.  Six 
teen  of  tin*  18  raw  data  files  can  be*  c'ditc'd  in  the  sc'iisc'  that  data  from  the-  raw 


Table  7.1:  Pre-Processor  Program  Main  Menu 


Module 

1  Rendezvous  Parameters  (TIGs) 

2  Rendezvous  Targeting  Parameters 

3  Sensor  Passes  (Time) 

4  Chaser  Vehicle  Halo  Orbit  Parameters 

5  Chaser  Vehicle  Initial  Conditions 

6  Chaser  Vehicle  T I  Error  Vectors 

7  Target  Vehicle  Halo  Orbit  Pai  ^meters 

8  Target  Vehicle  Initial  ('oiiditions 

9  Filter  Covariance  77  (’onditions 

10  Filter  Unmodeled  Acceleration  7 7  Conditions 

11  Filter  Measurement  Bias  and  Noi.se  77  Condilioiis 

12  Filter  Process  Noise  TI  Conditions 

13  Filter  Error  Plot  Scales  (Inertial) 

14  Filter  Error  Plot  Scales  (LVIR) 

15  Sensor  Biases 

16  Sensor  Errors 

17  Sensor  Limits 

18  Three- Body  System  Parameters 

19  Numerical  Integrator  Parameters 

20  Data  Storage  Flags 

21  Miscellaneous  Parameters 

99  Exit  Program 


1  KS 

data  file  will  be  stored  in  a  processed  data  file  at  a  selectable  interval.  In  this 
fashion  very  large  data  files  can  be  reduced  for  plotting  purposes.  These  raw 
data  files  are  never  written  over  or  deleted. 

The  data  file  containing  the  plotting  scales  for  the  estimation  error 
plots  can  be  completely  edited.  This  is  allowed  so  that  specific  areas  on  the  es¬ 
timation  error  plots  can  be  enlarged  without  re-executing  the  main  simulation 
or  editing  the  plotting  command  file.  The  last  raw  data  file  (for  miscellaneous 
parameter  storage)  can  not  be  edited.  Table  7.2  lists  the  main  menu  for  the 
post-processor  program  and  which  raw  data  files  are  available  for  editing.  I’lie 
type  of  data  stored  in  each  file  is  also  given. 
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Table  7.2:  Post-Processor  Program  Main  Menu 


Data  File 

Data  Type 

1 

chaser_env.dat 

Chaser  Environment 

2 

chaser_est.dat 

Estimated  States 

3 

chaser_guid.dat 

(^haser  Guidance  States 

4 

tgt_env.dat 

Target  Enviroiinieni 

5 

tgt_guid.dat 

'I'arget  Guidance  States 

6 

relmo.dat 

Relative  Motion  Data 

7 

errJnertial.dat 

Inertial  Estimation  Error 

8 

covJnertial.dat 

Inertial  Covariance 

9 

ccJnertial.dat 

Inertial  Correlation  Coefficient 

10 

errjvir.dat 

LVIR  Estimation  Error 

11 

covJvir.dat 

EVIR  Covariance 

12 

ccJvir.dat 

LVIR  Correlation  Coefficient 

13 

scales.dat 

Error  Plot  Scales 

14 

meas.dat 

MeasureiiK'iit  Data 

15 

residuals. dal 

Measurement  Residuals 

16 

updates. dat 

.Stale  I  pdates 

17 

foiii.dat 

Filter  Figiires-of  Merit 

99 

Exit  Program 

Chapter  8 


Conclusions 


A  guidance  law  has  been  developed  using  Ji-i  cont  rol  t  lieoi  \  wliic  li 
stabilizes  the  translunar  halo  orbit  in  the  circular  restricted  thrcx'  body  |)rob 
lem.  This  guidance  law  minimizes  the  position  deviation  liom  the  desiicul 
halo  orbit  plus  the  control  acceleration.  Linear  simulation  results  validated 
the  guidance  law.  Further,  the  halo  orbit  guidance  problem  has  been  formu¬ 
lated  in  the  frequency-domain.  Other  frequency- domain  design  techniques, 
such  as  Hoo  control  theory,  are  now  directly  applicable.  In  addition,  the  effect 
of  halo  orbit  frequency,  position  weighting  factor  and  the  amount  of  process 
or  measurement  noise  present  hcis  been  quantified. 

A  two  spacecraft  terminal  phase  rendezvous  targeting  law  has  been 
developed  which  is  valid  for  the  circular  restricted  three-body  problem.  'I'his 
targeting  law  was  demonstrated  using  a  small  radius  traiishmar  halo  orbil. 
In  addition,  several  characteristics  of  tlir<*(*-bo<ly  r»Mid»v,\ous  were  obseiAcil. 
First,  terminal  phase  rendezvous  in  the  circular  rest ii(  led  three  l)0(l\  |)iol> 
lem  is  three-dimensional:  out-of-plane  relativ«‘  motion  with  resix-i  t  to  the  lai 
get  vehicle’s  halo  orbit  occurs.  Further,  the  relative  trajectory  between  the 
chaser  and  target  vehicles  is  also  three-dimensional.  Secondly,  a  trade-off  ex¬ 
ists  between  targeting  law  error,  out-of-plane  relative  motion  and  the  total 
propulsion  cost  of  the  rendezvous  profile.  A  reiulezvous  profile  initiated  at 
the  top  or  bottom  of  the  halo  orbit  produces  less  targeting  law  error,  but  re¬ 
quires  greater  total  propulsion  cost  and  produces  larger  ont-of-plane  relative 
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motion;  a  rendezvous  profile  initiated  on  either  side  of  the  halo  orbit  requires 
less  total  propulsion  cost  and  produces  less  out-of-plane  relative  motion,  but 
has  greater  targeting  law  error.  Thirdly,  final  out-of-plane  erro’-s  generally 
dominate  both  downtrack  and  radial  errors.  Lastly,  a  minimum  total  propiil 
sion  cost  rendezvous  profile  exists  and  is  a  function  of  transfe:  lime,  iiu  line 
angle  and  initial  condition  angle. 

A  rendezvous  navigation  filter  capable  of  supplying  the  rendezvous 
targeting  law  with  chaser  vehicle  state  information  has  been  developed,  tis- 
timation  errors  converged  quickly  after  only  a  few  relative  range  and  line-of- 
sight  angle  measurements  were  incorporated.  The  placement  of  a  midcourse 
correction  maneuver  was  not  directly  dependent  on  navigation  errors  or  final 
position  error  sensitivities.  Hence,  the  inidcourse  correction  maneuvers  were 
placed  in  the  rendezvous  profile  based  on  time. 

F'our  total  cases,  two  rendezvous  profiles  for  the  translunar  planar  cir¬ 
cular  orbit  and  two  equivalent  rendezvous  profiles  for  the*  Hi  guided  translunar 
halo  orbit,  were  examined  to  demorrstrate  the  entire  guidance  and  na\  igatiun 
system  as  applied  to  the  two  spacecraft  terminal  phas(>  r('ud(‘/vous  scenario 
in  the  circular  restricted  three-body  |)robl<'m.  In  <'acli  <  a.se.  total  propulsion 
requirements  were  less  than  1..')  ni/s  and  final  miss  distance  was  less  than  iti 
m.  State  estimation  errors  were  small  and  showed  no  substantial  performance 
degradation  due  to  the  retidezvous  profile  geometry.  'I'lie  translunar  planar 
circular  orbit  displayed  the  same  characteristics  as  tim  H2  guided  translunar 
halo  orbit,  although  the  propulsion  requirements  and  final  miss  distance  were 
generally  larger  in  tln'TY^  guided  translunar  halo  orbit  demonstration. 

A  iiGilinear  simulation  package  for  terminal  pha.st'  rendezvous  guid¬ 
ance  and  navigation  system  evaluation  has  Ixvn  dev('lope<l.  The  pr«'-processor 
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and  post-processor  programs  are  completely  interactive;  the  main  simulation 
program  requires  no  interactive  input.  The  simulation  is  constructed  modu- 
larly  so  it  can  be  adapted  to  different  rendezvous  scenarios  in  any  three  l)o(ly 
system  with  a  minimum  of  modification.  The  .sinnilatioii  source  code  is  liillv 
documented  and  provides  processing  .status  messages  at  regular  intervals. 


Appendix  A 


State  Transition  Matrices 


A.l  Linearized  Equations  of  Motion 


Using  the  Ac  matrix  given  in  Eq.  (2.43),  the  state  transition  matrix, 


^g(^O),  can  be  computed  from  Eq.  (2.47)  where 


4>g(U0)  = 


^Gl  1  ^G\  2 

^G2l  ^G22 


(A.l) 


and 


*f*Gll 

‘PG\i\ 


<t>G\\i  4>aii2  0 
0  0 

0.667355 exp^  +0.6673.55  (*xp  • 


(A.2) 


-0.3;JI710cos(  1.862647/)  l.\.3| 

4>gi\2  =  -0.124821  cxp^''''^'^'+0.124S2Iexp-^''’'^‘"' 

+0.289317  sin(  1.862647/)  (A.4) 

4>gu*  =  -0.420595  exp''"'*^'^'  +0.420595  exp-*’'"'*^*""' 

+0.974880  sin(  1.862647/)  (A.5) 

4>Giis  =  0.0786672  exi)^''^'''  +0.0786672  exp-*'"*"'' 

+0.842666  co.s(  1 .862647/)  ( A. 6) 

<J>G\\9  =  cos(  1.786178/)  (.A  .7) 
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4>G\2 

4>GU\ 


4>G\22 


<i>Gl24 


<i>G\25 


<f>Gl29 


<t>GU\  <^G122  0 

<^G124  <^G125  0 

0  0  <i>G\29 

0.195181  exp^  -0.195181  exp--* 
+0.0844681  sin(  1.8626470 
0.123011  exp^  '“®^^'  +0.12301 1  exp"' 


(A.8) 

(.4.9) 


-0.246022  cos(  1 .862647/ )  (.4.10) 

-0.12301 1  exp^  -0.12301 1  cxp"^ 

+0.246022  cos(  1 .862647/)  (.4.11) 

-0.0775268  +0.0775268  exp-^*®*®^^' 

+0.716566  sin(  1.862647/)  (A. 12) 

0.559855  sin(  1.786 178/)  (A. 13) 


<^G21 

<i>G2U 

4>G2\2 

<t>G2\4 

<Ag215 


<i>G2\\  <t>G2\2  0 

<^G214  ^G215  0 

0  0  <Pg2\'J 

1 .440604  exp'^  '®^*^^'^'  - 1 .440604 

+0.623447  silt  (1.862617/) 

-0.269448 exp-*  -0.269448  cxp'^ 

+0. 538895  cos(  1.8626 17/) 

-0.907929exp2  ‘^’'’^‘ -0.907929 exp-^  '^*'''''' 

+ 1. 8 1 5857  cos(  1.86264  7/) 

0.169817exp^  '^^'' -0.169817 

-1. 569589  sin(  1.862647/) 

-1.786178sin(l. 786178/) 


(A.)  1) 


(.4.15) 


(.4.16) 


(.4.17) 


4>G2\^ 


(A.18) 

(.4.19) 
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^G22 

<i>G221 


<i>G22\  4>G222  0 

<i>G224  <f>G225  0 

0  0  <f>G229 

2.158677(  ‘  iX  /loi  999  «^rv“2.158677f 


0.421333 exp" -^0.42 1333 exp 
+0.157334  cos(l. 8626470 

<i>G222  =  0.265541  exp^  '^*®^'^' -0.265541  cxp-2  '^*®^^' 
+0.458252  sin(  1.8626470 
<f>G224  =  -0.265541  exp^  ’^*®^^' +0.265541  exp-" 

-0.458252  sin(  1.8626470 
4>G225  =  -0.167355exp'^ '***^”'-0.1673550X1)"" 

+ 1. 3347 10cos(  1.86264  70 
<i>G229  =  cos(  1.7861 780- 


(A.20) 

(A.21) 

(A.22) 

(A. 23) 

(A.24) 
(A. 25) 


A. 2  Rendezvous  Targeting  Law 

The  rendezvous  targeting  law  given  in  Eq.  (4.9)  is  dependent  on  the 
state  transition  matrix  computed  from  the  relative  state  equation,  Eq.  (4.6). 
Hence, 


<^11 

^12 

<t>u 

^21 

<f>22 

•  ^23 

^24 

A(; 

(A;  -  A,) 

0 

0 

•  <^3,3 

0 

A, 

0 

0 

^44  . 
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Pn 

V\2 

‘ei3 

<Pl4 

i*2\ 

022 

023 

024 

0 

0 

033 

034 

0 

0 

043 

044 

.(A.26) 


Solving  for  the  submatrices  of  interest  yields 


^111  <^112  0 

011  =  <^114  <^115  0  (A. 27) 

0  0  (^119 

<jiiii  =  0.667355  +0.667;J55exp-"''‘''‘'"'<'^-''> 

-0.334710 cos(l. 86264 7(<2  -  <i))  (A.28) 

4>u2  =  -0.124821  +0.124821  cxp"^ 

+0.289317sin(1.862647(<2  -  <i))  (A.29) 

<i>UA  =  -0.420595exp^'"‘^^^<'*-''’+0.420595exp-*'"“'^^<'»-‘') 

+0.974880 sin(1.862647(<2  -  <i))  (A.30) 

^,,5  =  0.078667exp^'®*^”<'^-‘'>  +0.078667  exp-^''‘“^'^’<‘^-'') 

+0.842666 cos(  1.8626  17(/2  -  /| ))  ( A.31 ) 

<;!>n9  =  cos(1.7S617S(/2  -  li))  (A. 32) 


4>l2l  4>li2  0 

012  =  *^124  <^125  0  (.\.33) 

0  0  <^>129 

<;i,2i  =  0.195181  exp^  -0.195181  exp“^  '''**‘^'''^-'‘ * 

+0.084468  sin(1.862647(/2  -  <i))  (A.34) 

<f>i22  =  0.123011  exp^'®**'"^'*"*'^ +0.12301  lexp~‘*'"'“'”<‘*“‘‘> 
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-0.246022  cos(1.862647(/2  -  <1))  (A.35) 

(f>i24  =  -0.12301  -0.123011 

+0.246022 cos(l. 86264 7(/2  -  ^i))  (A.36) 

(^125  =  -0.07  7527exp2'^®®^^<'^“''> +0.077527  exp--'’^*®'^*'^-''* 


+0.716566sin(1.862647(/2  -  <i))  (A.37) 

<^,29  =  0.559855sin(1.786178(/2  -  <i))  (A. 38) 

<^131  <^1,32  0 

013  *^134  <^135  0  (A. 39) 

0  0  <^139 

4>i3i  =  0.667335exp^‘'’'®^^^<'^"'‘>+0.667335exp~^’“®^^<'’“'’> 

-0.334710cos(1.862647(/2  -  /i))  -cos(f(/2  -  /,))  (A.40) 

<^132  =  -0.124820exp^'''^*®^^<‘^-‘'>  +0.124820(>xp-^'"^'’^'<'’-‘'> 

+0.289316 sin(  1.862647(/2  -  <1 ))  ( A.4 1 ) 

0134  =  -0.420595exp^  '^*‘'*'‘<'^-''>  +0.420595(>xp-^ 

+0.974880 siii(1.862647(/2  -  ^i))  (  A.  12) 

(^,35  =  0.078667('xp^''‘‘*^'‘‘*'^-'''+0.078667<‘xp-"-''''‘"‘'* 

+0.84  2666  (08(1. 86264  7(  <2  -  /i))  -  <o8(IV2  -  '1))  (  A.13) 

0139  =  cos(  1.786178(<2  - /i))  -  <os(f(/2  - /| ))  (A.  4  1) 


<^141  0142  0 

014  =  0144  014.S  0  (A. 45) 

0  0  0149 

0,4,  =  0.195181  fxp^  -0.195181  .•xp"''  > 

+0.084468 siii(1.862647(/2  -  /,))  -  -j:  siii(f{/2  -  /,))  (A.46) 
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0.12301 1  >  +0. 12301 1  > 

-0.246022 cos(l. 86264 7(/2  -  <i))  (.4.47) 

-0.123011  exp2'^*®^4«2-o)  _o.i2301  1  exp-^  '^*^""^-''> 
+0.246022 cos(1.862647(/2  -  <i))  (.\.18) 

-0.077526  ^0.077526 <‘xp“^  '''*^'"‘'^“''' 

+0.716565 sin(1.862647(<2  -  <i))  -  jsiii(f(<2  -  <i))  (A.49) 

0.559855 sin(1.786178(/2  -  <i))  -  |sin(f(<2  -  h))-  (A.50) 


Appendix  B 


Halo  Orbit  Guidance  Law  Data 


B.l  Representative  Guidance  Law 


The  mid-range  weighting  parameters  guidance  law  has  the  realization 


A-k 

Bk 

0 

where  Eqs.  (3.24),  (3.21)  and  (3.20)  yield,  respectively. 


-91.9710 

0.0011 

0 

0.0011 

-91.8698 
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-91.8680 

-3533.1400 

80.3422 

0 

76.7122 

-3518.3800 

0 

0 

0 

-1066.0300 

1.0000 

0 

0  ■ 

0 

1.0000 

0 

0 

0 

1.0000 

-14.1038 

44.3486 

0 

-2.0000 

0 

0 

0 

0 

-33.0702  . 

91.9740 

-0.0011 

0  ■ 

-0.0011 

91.8698 

0 

0 

0 

91.8680 

3525.7600 

76..5037 

0 

-76.7122 

.3516.1900 

0 

0 

0 

3516.0200  _ 

-1 1.7617 

|.5(>.8460 

0 

0 

0 

-.516.8190 

-1  1,1038 

12.3486 

0  ’ 

0 

0 

-33.0702 

(B.l) 


(B.2) 


(B.3) 


IB.I) 
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The  transfer  function  matrix  can  be  computed  from  the  state  space 
realization  using 


Ktf  =  Ck{sI  -  +  D„.  (B.5) 


Noting  Dfc  =  0,  Eq.  (B.5)  yields 


Ktf 


A’li  A'u  0 

0  0  A-23 


(B.h) 


where 


kn 

ki2 

^23 


-54333.24(5  +  1.06)(5^  +  97.625  +  3832.80) 
(52  +  98.925  +  3813.24  )(52  +  99.025  +  4573.70) 
162235.98(5  +  3.39)(52  +  91.285  +  3481.78) 
(52  +  98.925  +  3813.24)(.s2  +  99.025  +  4573.70) 
-166510.64(5+  11.49) 

(52  +  124.945  +  7104.12)' 


B.2  Representative  Time  Histories 

F'igure  B.l  shows  steady  state  time  historit's  for  the  syst(Mii  stat('s: 
Fig.  B.2  gives  the  estimation  error.  In  these  plots,  (a-e)  appl\  to  the  position 
states  while  (d-f)  apply  to  the  velocity  states.  Kigiin'  B.3  shows  ste<i(l\  state 
time  histories  for  the  coiit rol  acceleration  i'e<juired  to  |)iodim’  the  siea(l\  state 
positions  and  velocities. 


B.3  Representative  Halo  Orbits 

Figure  B.4  shows  the  front  view  of  halo  orbits  with  varying  amounts 
of  noise  included  in  the  simulation;  Fig.  B.5  giv<>s  the  side  view  for  the  sane 


cases. 
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B.4  Noise  Statistics 

Figure  B.6  shows  the  average  standard  deviation  of  tiu'  position  and 
velocity  states  as  a  function  of  the  standard  deviation  of  the  process  noise.  In 
both  figures,  the  fixed  power  spectrum  of  the  measurement  noise  is  assumed 
to  be  2.665  m*/s  in  each  axis. 

B.5  Halo  Orbit  Frequency  Parametric  Study 

Figures  B.7  and  B.8  show  the  geometry  and  orientation  of  the  halo 
orbit  as  a  function  of  non-dimensional  frequency.  Figures  B.9  and  B.IO  give 
the  velocity  required  as  a  function  of  non-dimensional  frequency  in  the  ^  and 
C-axes,  respectively;  Fig.  B.  1 1  applies  to  fii<‘  total  velocity  r(‘()iiirf'd.  In  all 
these  figures,  (a-c)  apply  to  cases  with  constant  wc'iglitiiig  in  llic  //  ax is  and 
(d-f)  apply  to  cases  with  constant  WfMgliting  in  the  a.\is. 


B.6  Control  Weighting  Factors  Parametric  Study 

F'igure  B.r2  shows  the  effect  of  the  weighting  factor  in  the  r/-axis 
for  constant  values  of  halo  orbit  frequency  and  weighting  factor  in  the  (,'-axis. 
Figure  B.13  shows  the  effect  of  the  weighting  factor  in  the  (-axis  for  constant 
values  of  halo  orbit  frerjiiencv  and  weig'.iting  factor  in  the  //-axis. 
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Figure  B.4:  Halo  Orbit  (leoiuetry  a.s  a  Function  of  Nois< 
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Figure  B.8;  Halo  Orbit  Orientation  as  a  Function  of  Frequency 
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Figure  B.!):  VVlocity  in  llif^  Axis 


Frequency  (Non'<limeneiontl) 


Frequency  (Non-dimoitional) 


Figure  IMl:  'I’otal  V<‘lority  HiHiuired 


Appendix  C 


Rendezvous  Navigation  Filter  Monte  Carlo  Data 

C.l  Individual  Error  Vector  Plots 

Figures  C.l  -  C.13  show  the  estimation  error  (solid  line)  from  the 
rendezvous  navigation  filter  of  the  chaser  vehicle’s  inertial  state  for  each  of 
the  13  errors  vectors  selected.  Range  and  linc-of-sight  angle  measurements 
were  incorporated  into  the  rendezvous  navigation  filter  once  every  minute.  In 
these  plots,  (a-c)  apply  to  the  position  estimation  error  while  (d-f)  apply  to 
the  velocity  estimation  error.  The  1<t  standard  deviation  of  the  state  (dotted 
lines)  are  also  included  on  the  plots. 

C.2  Sensitivity  Matrix  Data 

Figures  (’.14  -  (Mfi  show  the  time  liislori('s  of  the  individual  eleiiients 
of  the  sensitivity  matrix.  These  figures  apply  to  clianges  in  (h)wntra(  k.  ont- 
of-plane  and  radial  velocity,  re'^pectively. 


Fifiure  C.2:  Chaser  V^ehicle  State  Kstimation  Krror  (Krror  Vector  2) 


Figuro  C.3:  Chaser  Vehicle  State  Kstiination  Krror  (F.rror  Vector  3) 


Figure  C.8:  Chaser  Vehicle  Stat«*  Estimation  Error  (Error  Vet  tor  S) 


(b)  YL  Poi  EiTor  (m)  vs.  Time  (Kts) 


Figure  C.l.J:  Chaser  Vehicle  J; 


t*  Fsliiiiatiori  Krror  (Krror  Vector  13) 


Appendix  D 


Applications  Monte  Carlo  Data 

D.l  Planar  Circular  Orbit  (Profile  1) 

Table  D.l  summarizes  the  propulsion  requirements  for  the  individual 
simulations  of  the  planar  circular  orbit  for  profile  1.  Table  D.2  gives  the  final 
miss  distance  data  for  these  same  simulations. 


D.2  Planar  Circular  Orbit  (Profile  2) 

Table  I). 3  summarizes  the  propulsion  requireiiienis  for  the  indiv  idual 
simulations  of  the  planar  circular  orbit  for  profile  2.  Table  1).  I  i>i\<‘s  the  liiial 
miss  distance  data  for  these  same  simulations. 


D.3  7(2  Guided  Halo  Orbit  (Profile  1) 

Table  D.5  summarizes  the  propulsion  requirements  for  the  individual 
simulations  of  the  guided  halo  orbit  for  profile  1.  Table  1).6  gives  the  final 
miss  distance  data  for  these  same  simulations. 


D.4  7i2  Guided  Halo  Orbit  (Profile  2) 

Table  D.7  summarizes  the  propulsion  requirements  for  the  individual 
simulations  of  the  H2  guided  halo  orbit  for  profile  2.  Table  D.S  gives  tin-  final 
miss  distance  data  for  these  same  simulations. 


I(i2 


Table  D.l:  Planar  Circular  Orbit  Propulsion  Requirement  (IC=46.08  deg, 
INC=22.85  deg,  TT=9  hrs) 


Error  Rendezvous  Maneuver  (m/s)  R(’S 


Vector 

TI 

MCCl 

MCC2 

Mcc:\ 

MCCA 

TF 

'I'otal 

Fuel  (lbs) 

Reference 

0.686 

0.000 

0000 

0.000 

0.000 

0  491 

i .  177 

5  776 

1 

0.686 

0  032 

0.009 

0.016 

0.018 

0.487 

1.247 

6.119 

2 

0.686 

0.016 

0003 

0.014 

0.017 

0,500 

1  236 

6.065 

3 

0.686 

0.055 

0.006 

0.010 

0.017 

0.454 

1.228 

6.026 

4 

0.686 

0.046 

0.008 

0.009 

0.005 

0.453 

1.207 

5.923 

5 

0.686 

0.036 

0.013 

0.013 

0.021 

0.504 

1.272 

6.242 

6 

0.686 

0.066 

0.014 

0.012 

0.010 

0.523 

1  310 

6.428 

7 

0.686 

0.033 

0.018 

0.010 

0.006 

0.477 

1.230 

6.036 

8 

0.686 

0.013 

0.010 

0.014 

0.022 

0.478 

1.223 

6.001 

9 

0.686 

0.039 

0.002 

0.007 

0.016 

0.492 

1.242 

6.094 

10 

0.686 

0.056 

0.018 

0.005 

0.014 

0.490 

1.268 

6.222 

11 

0.686 

0.047 

0.009 

0.007 

0.008 

0.496 

1.252 

6.143 

12 

0.686 

0.072 

0.006 

0.010 

0.011 

0.554 

1.339 

6.570 

13 

0.686 

0.018 

0.012 

0.007 

0.019 

0.504 

1.246 

6.114 

Average 

Magnitude 

0.686 

0.041 

0.010 

0.010 

0.014 

0.493 

1.254 

6.153 
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Table  D.2:  Planar  Circular  Orbit  Final  Miss  Distance  (IC=46.08  deg 
INC=22.85  deg,  TT=9  hrs) 


Error  Final  Miss  Distance  (m) 


Vector 

Downtrack 

Out-of- Plane 

Radial 

Total 

Reference 

0.000 

-0.001 

0.000 

0.001 

1 

2.615 

14.744 

-3.113 

15.295 

2 

-12.508 

-2.864 

-9.094 

15.728 

3 

-3.617 

-13.090 

-6.687 

15.138 

4 

-5.347 

2.661 

-12.409 

13.771 

5 

9.999 

-5.066 

8.7.39 

14.213 

6 

2.108 

-1.856 

-5.405 

7.-566 

7 

-3.099 

-6.349 

-13.804 

15.-507 

8 

-0.7.30 

9.458 

9.759 

13,610 

9 

-0.977 

-7.665 

- 1,359 

8.N72 

10 

10.051 

-0.343 

-2.66.3 

10,103 

11 

4.840 

5.942 

8.-54  0 

11,171 

12 

-3.-58 1 

9.121 

3.776 

10.-502 

13 

4.469 

11.498 

2.437 

12.574 

Average 

Magnitude 

4.919 

7.204 

6.983 

12.666 
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Table  D.3: 
INC=22.85 

Planar  Circular  Orbit  Propulsion  Requirements  (IC= 
deg,  TT=9  hrs) 

180.00  deg. 

Error 

Rendezvous  Maneuver  (m/s) 

RCS 

Vector 

TI 

MCC\ 

MCC2 

Mca 

MCC4 

TF 

Total 

Fuel  (lbs) 

Reference 

0.680 

0.000 

0.000 

0.000 

0.000 

0.643 

1.322 

6,487 

1 

0.680 

0.035 

0.005 

0.008 

0.010 

0.610 

1.348 

6.614 

2 

0.680 

0.009 

0.008 

0.010 

0.017 

0.649 

1.372 

6.732 

3 

0.680 

0.058 

0.013 

0.007 

0.013 

0.665 

1.435 

7  041 

4 

0.680 

0.042 

0.011 

0.009 

0.020 

0.662 

1.423 

6.982 

5 

0.680 

0.036 

0.012 

0.004 

0.009 

0.640 

1.381 

6.776 

6 

0.680 

0.060 

0.007 

0.012 

0.023 

0.634 

1  115 

6  913 

7 

0.680 

0.039 

0.009 

0.005 

0.013 

0.622 

1.368 

6.712 

8 

0.680 

0.014 

0.006 

0.003 

0.012 

0.632 

1.346 

6.605 

9 

0.680 

0.038 

0.010 

0.011 

0.006 

0.628 

1.372 

6.732 

10 

0.680 

0.059 

0.004 

0.007 

0.020 

0.650 

1.420 

6.968 

11 

0.680 

0.049 

0.009 

0.008 

0.013 

0.611 

1.370 

6.722 

12 

0.680 

0.073 

0.010 

0.005 

0.012 

0.619 

1.399 

6.865 

13 

0.680 

0.019 

0.007 

0.015 

0.014 

0.616 

1.350 

6.624 

Average 

Magnitude 

0.680 

0.041 

0.009 

0.008 

0.014 

0.634 

1.384 

6.791 
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Table  D.4:  Pltinar  Circular  Orbit  Final  Miss  Distance  (1('=  180.00  deg, 
INC=  22.85  deg,  TT=9  hrs) 


Error 

Final  Miss  Distance  (m) 

Vector 

Downtrack  Out-of-Plane 

Radial 

Total 

Reference 

0.000 

0.000 

-0.001 

0.001 

1 

6.641 

12.333 

-2.685 

14.262 

2 

-14.683 

-8.169 

-3.234 

17.111 

3 

4.325 

-4.754 

8.704 

10.820 

4 

-6.652 

-3.291 

-7.559 

10.593 

5 

1.331 

-12.208 

7.675 

14.481 

6 

13.045 

14.898 

-3.024 

20.032 

7 

-12.729 

-15.405 

-1.809 

20.065 

8 

4.467 

2.851 

6.503 

8.389 

9 

-21.063 

-15.329 

-4.528 

26.441 

10 

6.461 

-5.222 

-5.814 

10.140 

11 

-5.767 

-1.066 

-9.905 

11.511 

12 

-1.938 

6.628 

6.501 

9.485 

13 

6.021 

24.977 

5.973 

26.378 

Average 

Magnitude 

8.08: 

9.779 

5.686 

15.362 

Table  D.5:  W2  Guided  Halo  Orbit  Propulsion  Requirements  (IC=46.08  deg, 
TT=9  hrs) 


7.071 
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Table  D.6:  Guided  Halo  Orbit  Final  Miss  Distance  (1( '=46.08  deg.  TT=9  » 

hrs) 


Error 

Vector 

Final  Miss  Distance  (m) 

Downtrack 

Out-of- Plane 

Radial 

Total 

Reference 

-7.673 

19.394 

0.124 

20.857 

1 

-10.865 

26.681 

3.703 

29.046 

2 

-15.263 

12.383 

22.215 

29.662 

3 

-11.777 

10.377 

-25.045 

29.558 

4 

-7.798 

4.792 

17.163 

19.451 

5 

-2.458 

9.457 

-3.898 

10.520 

6 

4.001 

32.472 

-4.514 

33.027 

7 

-10.714 

20.000 

-9.875 

24.745 

8 

-9.348 

8.864 

4.978 

13.811 

9 

-5.351 

24.825 

-4.748 

25.835 

10 

-21.348 

31.346 

-1.280 

37.947 

11 

-5.276 

36.849 

8.084 

38.093 

12 

-1.354 

24.097 

10.064 

26.149 

13 

-13.530 

12.465 

-5.375 

19.166 

Average 

Magnitude 

9.160 

19.585 

9.303  25.9‘24 

I 
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Table  D.7:  'H2  Guided  Halo  Orbit  Propulsion  Requirements  (1C’=  180.00  deg, 
TT=9  hrs) 


Error 

Rendezvous  Maneuver  (m/s) 

RCS 

Vector 

TI 

MruT 

MCC2 

A/CC3 

MCC4 

TF 

Total 

Fuel  (lbs) 

Reference 

0.704 

0.008 

0.008 

0.009 

0.013 

0.672 

1.413 

6.933 

1 

0,704 

0.029 

0.020 

0.003 

0.032 

0.655 

1.442 

7.076 

2 

0.704 

0.029 

0.007 

0.012 

0.037 

0.687 

1.476 

7.242 

3 

0.704 

0.055 

0.012 

0.017 

0.026 

0.694 

1.509 

7.404 

4 

0.704 

0.048 

0.009 

0.024 

0.003 

0.691 

1.481 

7.267 

5 

0.704 

0.051 

0.011 

0.014 

0.014 

0.664 

1.457 

7.149 

6 

0.704 

0.059 

0.013 

0.021 

0.016 

0.651 

1.465 

7.188 

7 

0.704 

0.025 

0.009 

0.010 

0.026 

0.650 

1.424 

6.987 

8 

0.704 

0.007 

0.016 

0.011 

0.013 

0.656 

1.408 

6.909 

9 

0.704 

0.040 

0.021 

0.012 

0.021 

0.663 

1.460 

7.164 

10 

0.704 

0.064 

0.013 

0.010 

0.024 

0.684 

1.498 

7.350 

11 

0.704 

0.054 

0.012 

0.009 

0.013 

0.640 

1.433 

7.031 

12 

0.704 

0.070 

0.017 

0.007 

0.029 

0.640 

1.468 

7.203 

13 

0.704 

0.013 

0.015 

0.010 

0.013 

0.646 

1.400 

6.870 

Average 

Magnitude 

0.704 

0.042 

0.013 

0.012 

0.021 

0.663 

1.455 

7.139 
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Table  D.8:  H2  Guided  Halo  Orbit  Final  Miss  Distance  (1(’=  180.00  deg,  TT=9 
hrs) 


Error 

Vector 

Final  Miss  Distance  (m) 

Downtrack 

Out-of- Plane 

Radial 

Total 

Reference 

-4.834 

-3.191 

-0.885 

5.859 

1 

-2.708 

-5.606 

0.918 

6.293 

2 

4.658 

1.088 

-2.222 

5.274 

3 

-4.074 

-11.543 

0.630 

12.257 

4 

5.344 

3.019 

-4.398 

7.551 

5 

-7.830 

-18.871 

3.705 

20.764 

6 

-4.654 

2.036 

7.515 

9.070 

7 

-2.251 

8.457 

0.655 

8.776 

8 

-12.955 

-10.494 

-1.436 

16.734 

9 

6.421 

-26.408 

19.593 

33.504 

10 

-13.834 

-16.410 

5.362 

22.122 

11 

-1.925 

-11.195 

-14.518 

18.434 

12 

-5.174 

-3.076 

3.882 

7.162 

13 

-4.848 

-14.685 

0.6.39 

15.478 

Average 

Magnitude 

5.898 

10.222 

5.036 

1 1.109 
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